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*0 ' Abstract 

Given a time symmetric initial data set for the vacuum Einstein field equations which 
is conformally flat near infinity, it is shown that the solutions to the regular finite initial 
value problem at spatial infinity extend smoothly through the critical sets where null infinity 
touches spatial infinity if and only if the initial data coincides with Schwarzschild data near 



?H , infinity. 
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^ ■ 1 Introduction 

i ! _ 

[ — . A vacuum spacetime (A4,g^ u ) is said to be asymptotically simple — cfr. [161 1191 [25] — if there 

exists a smooth, oriented, time-oriented, causal spacetime (M, g^ v ) and a smooth function H (the 
\Q . conformal factor) on M. such that: 

O ' 

' (i) A4 is a manifold with boundary J r , 



o 



(ii) S > on M \ J and S = 0, dS ^ on J, 

(iii) there exists a conformal embedding $ : Ai — ► M. \ J? such that 



(iv) each null geodesic of {ftA,g^) acquires two distinct endpoints on J> . 

In this definition and throughout this article the word smooth is used as a synonym for C°° . The 
point (iv) in the definition of asymptotic simplicity excludes black hole spacetimes. In order to 
incorporate this class of spacetimes into the framework, the definition has to be modified. A way 
of doing this is through the notion of weak asymptotic simplicity — cfr. |25l I27j . This point will 
not be of relevance in our considerations. 

The definition of asymptotic simplicity was introduced in [24j as a way of characterising the 
asymptotic behaviour of spacetimes describing isolated systems. As such, it allows to reformulate 
questions concerning the asymptotic decay of physical fields in terms of its differentiability at the 
conformal boundary. In particular, asymptotic simplicity implies the peeling behaviour on the 
Weyl tensor of the physical spacetime (A4,g^) — see e.g. |25) . 

It has been long debated whether the condition on the smoothness of the conformal boundary 
is a too strong requirement to be imposed on a spacetime describing interesting physical phenom- 
ena. In [SJ it has been shown that there exists a large class of non-trivial asymptotically simple 
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spacetimes. These spacetimes have the peculiarity of being isometric to the Schwarzschild space- 
time in a neighbourhood of spatial infinity. Given this particular state of affairs, the question is 
now: is this the only possible way to construct non-trivial asymptotically simple spacetimes? 

The Cauchy problem in General Relativity provides a systematic approach to address the 
existence of more general classes of asymptotically simple spacetimes. In the spirit of asymptotic 
simplicity, it is natural to work with objects and equations which are defined in the conformally 
rescaled spacetime {M.,g^ v ). The conformal Einstein field equations introduced in [9j[8l[T0] and 
extensions thereof — see [TH [T51 HH HH] — provide a suitable system of equations and unknowns. 
These conformal field equations are formally regular at the points where the conformal factor 
S vanishes, and such that a solution of them implies a solution to the vacuum Einstein field 
equations. In particular, the conformal Einstein field equations have been used to prove a semi- 
global existence and stability result for hyperboloidal data [12] — see also [21]. The construction 
of asymptotically simple spacetimes in [5] makes use of this result. 

Cauchy data for asymptotically simple spacetimes is prescribed on hypersurfaces which are 
asymptotically Euclidean. The compactification of an asymptotically Euclidean hypersurface 
includes a singled out point for each asymptotic end — the corresponding "point at infinity". 
These points become the spatial infinities of the development of the Cauchy data. Note however, 
that the definition of asymptotic simplicity makes no reference to these points; the reason being 
that if the spacetime has non-zero ADM mass then spatial infinity is a singular point of the 
conformal geometry — see e.g. [13] . 

It has long been acknowledged that the main obstacle in the construction of asymptotically 
simple spacetimes from Cauchy initial data is the lack of a detailed understanding of the structure 
of the Einstein field equations at spatial infinity — cfr. [TTJ Q2J . The difficulties in the analysis 
of the structure of spatial infinity arise from the singular nature of this point with respect to the 
conformal geometry — as mentioned in the previous paragraph. In particular, the Weyl tensor 
is singular at spatial infinity if the ADM mass is non-zero. Thus, one needs to obtain a better 
representation of spatial infinity. 

The regular finite initial value problem near spatial infinity introduced in [15j provides a 
powerful tool for analysing the properties of the gravitational field in the regions of spacetime 
"close" to both spatial and null infinity. This initial value problem makes use of the so-called 
extended conformal Einstein field equations and general properties of conformal structures. It is 
such that both the equations and the data are regular at the conformal boundary — the regular 
finite initial value problem at spatial infinity. 

Whereas the standard compactification of spacetimes put forward in the definition of asymp- 
totically simple spacetimes considers spatial infinity as a point, the approach used in [15] repre- 
sents spatial infinity as an extended set with the topology [—1,1] x § 2 . This so-called cylinder 
at spatial infinity is obtained as follows: starting from an asymptotically Euclidean initial data 
set for the Einstein vacuum equations, one performs a "standard" conformal compactification to 
obtain a compact manifold, S, with a singled out point, i, representing the infinity of the initial 
hypersurface — for simplicity it is assumed there is only one asymptotic end. In a second stage, 
the point i is blown-up to a 2-sphere. This blowing-up of i is achieved through the lifting of a 
neighbourhood of i to the bundle of space-spinor dyads. In the final step, one uses a congruence 
of timelike conformal geodesies to obtain Gaussian coordinates in a neighbourhood of the initial 
hypersurface. The conformal geodesies are conformal invariants that can be used to construct 
a gauge which renders a canonical class of conformal factors for the development of the initial 
data. These conformal factors can be written entirely in terms of initial data quantities. Hence, 
the location of the conformal boundary is known a priori. The conformal boundary described 
by the canonical conformal factor contains a null infinity with the usual structure and a spatial 
infinity which extends in the time dimension — so that one can speak of the cylinder at spatial 
infinity. The sets {±1} x S 2 will be called critical sets as they can be regarded as the collection of 
points where null infinity "touches" spatial infinity. Null infinity and spatial infinity do not meet 
tangentially at the critical points. As a consequence, part of the propagation equations implied 
by the conformal field equations degenerates at these sets. The analysis in [15] has shown that, as 
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a result, the solutions to the conformal field equations develop certain types of logarithmic singu- 
larities at the critical sets. These singularities form an intrinsic part of the conformal structure. 
It is to be expected that the hyperbolic nature of the propagation equations will propagate these 
logarithmic singularities. Consequently, they will have an effect on the regularity of null infinity. 
However, there is not proof of this yet. 

The analysis carried out in [15] — restricted to the developments of time symmetric data 
with an analytic conformal compactification at infinity — rendered an infinite set of necessary 
conditions on the initial data for the spacetime to extend analytically through the critical sets. 
These conditions can be formulated in terms of the Cotton tensor of the initial metric and its 
higher order derivatives. If these conditions are not satisfied at some order then the solutions to 
the conformal Einstein equations develop singularities of a very definite type at the critical sets. 
It is important to point out that these singularities are associated to structural properties of the 
principal part of the evolution equations. 

The ideas and techniques of [15 can be implemented in a computer algebra system. In |30j this 
approach has been used to analyse a certain type of asymptotic expansions of the development of 
initial data sets which are conformally flat in a neighbourhood of infinity. The rationale behind the 
use of conformally flat data sets is that they satisfy the regularity conditions of [TS] automatically. 
They constitute the simplest type of initial data sets on which the methods of |15j can be applied 
to obtain non-trivial results. Although analytically simple, conformally flat initial data sets give 
rise to spacetimes of full complexity, and indeed, particular examples are used routinely in the 
simulation of head-on black hole collisions — see e.g. [U [4j [22] . The key finding of [30] was to 
show the existence of a further type of logarithmic singularities at the critical sets. This new 
class of obstructions to the smoothness of null infinity arises from the interaction of the principal 
part (which is singular at the critical sets) and the lower order terms in the conformal Einstein 
field equations. An important observation arising from the analysis in [30j is that particular 
subsets of these singularities can be eliminated by setting certain pieces of the initial data to 
zero. Proceeding in this way it is possible to gain insight into the algebraic structure of the 
conformal field equations at spatial infinity, and a very definite pattern is observed. One has the 
following conjecture. 

Conjecture. If an initial data set for the Einstein vacuum equations which is time symmetric and 
conformally fiat in a neighbourhood of infinity yields a development with a smooth null infinity, 
then the initial data set is exactly Schwarzschildean in the neighbourhood. 

In this article we make definite progress towards a proof of this conjecture. The analysis 
presented in this article is concerned with the behaviour of the solutions of the conformal Einstein 
field equations at the critical sets. The main result of this article is the following: 

Theorem. Consider a time symmetric initial data set for the Einstein vacuum field equations 
which is conformally flat near infinity. The solution to the regular finite initial value problem at 
spatial infinity is smooth through the critical sets if and only if the data is exactly Schwarzschildean 
in a neighbourhood of infinity. 

The result presented here falls short of providing a proof of the conjecture for it is a priori not 
clear that smoothness at null infinity follows from the smoothness at the critical sets — although 
the expectation is that this will be the case. Conversely, the connection between a singular 
behaviour at the critical sets and non-smoothness at null infinity is, as yet, not fully understood. 
It is expected that estimates of the type constructed in [18] should provide a way of linking the 
asymptotic expansions obtained by the methods of [15j with actual solutions to the conformal 
field equations. A generalisation of these estimates to the setting of a linear field propagating on 
a curved background has been constructed in [35] . 

In order to bring our main result into a wider context, it is pointed out that the Schwarzschild 
spacetime is the only static spacetime admitting conformally flat slices — see e.g. [3] [19]. Thus, 
our main result together with the analysis in [29j of asymptotic expansions for more general 
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classes of time symmetric data suggests that it should be possible to prove a generalisation in 
which smoothness at the critical sets implies staticity in a neighbourhood of spatial infinity. The 
situation for initial data sets with a non-vanishing second fundamental form is not as clear cut 
as one loses the analyticity of the conformal compactification — cfr. [6]. In any case, initial data 
sets which are stationary in a neighbourhood of infinity are expected to play a privileged role 
—see e.g. [3H [32j [33] . 

The analysis leading to the main theorem provides very detailed information about the be- 
haviour of the solutions to the conformal Einstein field at spatial infinity. In particular, it will be 
shown that these singularities appear at higher orders than what it is to be expected from general 
arguments. Indeed, there are some structural properties of the conformal field equations (some 
cancellations) which make the solutions more regular than what they are expected to be. This 
explains why the first analysis carried out in [20] did not observe any singularities, and one had to 
undertake the higher order expansions of 30J. It is certainly desirable to obtain a deeper under- 
standing of the structures responsible for this behaviour. It is important to emphasise that these 
cancellations make the analysis much more computationally challenging than what it otherwise 
would be. 

Outline of the article 

The present work is based on the analysis of [15] — cfr. also [19] . It makes use of a number 
of structures and constructions which are not very standard. Therefore, in order to make our 
discussion accessible, it will be necessary to introduce a certain amount of notation and definitions. 
This should not be taken as a comprehensive overview. The reader is, in any case, referred to 
reference [15] for complete details. I have striven, as much as it is possible, to respect the original 
conventions of [15 . The main difference is the use of capital Latin letters to denote spinorial 
indices. 

Section [2] introduces some basic notational conventions and presents the notion of asymptot- 
ically Euclidean data in the form that will be used here. 

Section [3] provides a brief overview of time symmetric, conformally flat initial data sets and 
introduces the notion of data which is Schwarzschildean up to a certain order. 

Section 0] provides a summary of the blowing-up of the point at spatial infinity. The blow-up 
is realised by the introduction of a certain manifold C a which is a subset of the bundle of spin- 
frames. It is also shown how one can introduce a calculus on this manifold, and the form that 
normal expansions acquire when lifted up to this manifold. The contents of this section follow 
reference [15 a , of which they constitute a very terse summary. 

Section [5] discusses the extended conformal field equations and the evolution equations they 
imply when written in a conformal Gaussian gauge system — the F-gauge. It discusses some 
structural properties of these evolution equations and introduces the notion of transport equations 
at the cylinder at spatial infinity. 

Section [6] discusses the Schwarzschild spacetime in the F-gauge. This discussion will be a 
constant reference point. 

Section presents the decomposition of the linear transport equations into spherical harmon- 
ics. It introduces the notion of a spherical harmonic sector and discusses the discrete symmetries 
implied by the time reflection symmetry of the spacetime. It also presents a procedure to solve 
the hierarchy of transport equations. This procedure is different to the one given in |15j as it 
exploits a certain analogy with the Maxwell equations — cfr. [34] , 

Section [5] presents a systematic analysis of the properties of the solutions of the transport 
equations for initial data sets which are Schwarzschildean up to a certain order. It briefly reviews 
the computer algebra calculations of reference [30 . It is shown that the solutions to these trans- 
port equations are more regular than what one would a priori expect, but that neverthelesss, 
they eventually develop certain logarithmic singularities at the critical sets. The order at which 
these singularities arise is determined. The key results of this section are based on very lengthy 
computer algebra calculations. Only qualitative aspects of these calculations are presented. 
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In section|9l the results of the previous sections are summarised in theorem[2] This theorem is 
a more precise version of the main result discussed in the introductory paragraphs. Some further 
discussion is given, in particular regarding the potential generalisation to initial data sets which 
are not conformally flat. 



2 Notation and conventions 

This article is concerned with the asymptotic properties of spacetimes {M^^u) solving the Ein- 
stein vacuum field equations 

R„ v = 0. (l) 

The metric g^ v will be assumed to have signature (+, — , — , — ) and fi, v, . . . are spacetime indices 
taking the values 0, . . . , 3 . The spacetime (M., g^) will be thought of as the development of some 
time symmetric initial data prescribed on an asymptotically Euclidean Cauchy hypersurface S. 
The time symmetric data on S are given in terms of a 3-metric hij of signature (— , — , — ). Due to 
the requirement of time symmetry one has that the other piece of the data, the second fundamental 
form of iS vanishes: Xij = 0- The letters i, j, k, . . . will be used as spatial tensorial indices taking 
the values 1, 2, 3. In this case the Einstein vacuum field equations imply the constraint equation 

r = 

where f denotes the Ricci scalar of the metric hij. For simplicity, only one asymptotically flat 
end will be assumed. In the asymptotically flat end it will be assumed that coordinates can 
be introduced such that 

hij = - (l + y^j Sij + O (Jj^j , as \y\ oo, 

with \y\ 2 — (y 1 ) 2 + (y 2 ) 2 + (y 3 ) 2 and m a constant — the ADM mass of S. In addition to these 
asymptotic flatness requirements, it will be assumed that there is a 3-dimensional, orientable, 
smooth compact manifold (S,hij), a point i € S, a diffcomorphism $ : S \ {i} — ► S and a 
function fi G C 2 (S) nC°°(S\ {i}) with the properties 

n(i) = 0, DjCl(i)=0, DjD k n(i) = -2h Jk (i), (2a) 
n > on S \ {i}, (2b) 
hij = Vl 2 §Jiij. (2c) 

The last condition shall be, sloppily, written as hij = Vl 2 hij — that is, S \ {i} will be identified 
with S. Under these assumptions (S, hij) will be said to be asymptotically Euclidean and regular. 
Suitable punctured neighbourhoods of the point i will be mapped into the asymptotic end of 5. 
It should be clear from the context whether i denotes a point or a tensorial index. 



3 Conformally flat, time symmetric initial data sets 

All throughout this article it will be assumed that one has initial data sets which are conformally 
flat in a neighbourhood B a C S, a > 0, of i G S. The expressions "near infinity" and "in a 
neighbourhood of infinity" used in the introductory paragraphs should be understood in this 
sense. 

The Hamiltonian constraint together with the boundary conditions ([2a j) -([2c |) imply on B a (i) 
the Yamabe equation 

Ai? = -47n5(i), tf = fT 2 , (3) 

where S(i) denotes the Dirac delta distribution with support on i, and A is the flat Laplacian. 
Let x l denote Cartesian coordinates on B a such that x l (i) — 0. On B a \i the physical 3-metric 
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hij is given by hij = — i? 5^, where <£y is the standard Euclidean metric in Cartesian coordinates. 
If a > is suitably small, the solutions of can be written on B a in the form 

*=i + W, (4) 



where \x\ = v/ (a; 1 ) 2 + (a; 2 ) 2 + (a; 3 ) 2 . The function W satisfies 

AW = 0, (5) 

so that W is a harmonic and analytic function. From general considerations on asymptotics, one 
has that 

where m is the ADM mass of the time symmetric initial data set (S, hij) — cfr [T21 HS1 US] ■ More 
generally, one has that 

oo 

W=j + ^n-* k x n ---x ik , (6) 

with w^.-.if. denoting constant tensors which are totally symmetric and (5-tracefree so that the 
expression Wi 1 ...i k x tl ■■■x lk is a homogeneous harmonic polynomial of degree k. Conversely, 
due to the analyticity of the solutions to the Laplace equation |5j), given a sequence {ttj^...^}, 
k = 0, . . . , oo, of constant totally symmetric and trace-free constant tensors (the germ of W) such 
that the sum 

oo 

y t w ii---i k x 1 ■ ■ • X k , 
fc=l 

converges in B a for a given a > 0, then the sequence defines a unique solution of the Yamabe 
equation ((3]). 

The following well-known observation will simplify our analysis. Consider the inversion y % = 
x l /\x\ 2 mapping B a to the asymptotic end R 3 \ Bi/ a (0), and let <fi = \x\d. If 

$ = i — r + -jr + w i x i + w tJ x l x j + 0(\x\ 3 ), 
\x\ I 

then it follows readily that 

11 m I i I i i ^ ( I 

%l M 3 M 5 V 1 2/ 1 

Now, consider the translation z J = j/* + c*, with c l a constant vector in R 3 . A direct calculation 
shows that 

1 If c,z l „ / 1 \\ 1 1 / „ / 1 



Ivl 1*1 V 1 M 2+0 Ul 2 JJ' M 3 kl 3 l 1 + °VNI 



from where it follows that 

777 



Thus, the dipolar term can be removed by setting Cj = Wi/m — that is, by choosing appropriately 
the "centre of mass". Finally, an inversion of the form z % = x l /\x\ 2 , shows that one can restrict 
the analysis to solutions to equation (J3j> which in B a {i) take the form 

^=^- + ^+W ij xV+0(\x\ 3 ). 

In what follows it will always be assumed that 7? is given in the above form. In the above 
representation, data which is exactly Schwarzschildean on B a is characterised by 

7/7 
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More generally, we have the following. 



(i 



Definition. A time symmetric initial data set which is conformally flat in a suitable neighbour- 
hood B a of infinity will be said to be Schwarzschildean up to order p, if and only if 

^ = R + ? + ^•••w^ 1 • • • xlp ' +1 + o(\x\ p - +2 ), 

that is, if 

W -m/2 = <D(\x\ p ' +1 ). 



4 The Manifold C a 

In [15] a representation of the region of spacetime close to null infinity and spatial infinity has 
been introduced — see also the comprehensive discussion in [T5]. The standard representation 
of this region of spacetime depicts i° as a point. In contrast, the representation introduced in 
1 1 5 j depicts spatial infinity as a cylinder — the cylinder at spatial infinity. This construction is 
briefly reviewed for the case of time symmetric initial data sets which are conformally flat in a 
neighbourhood B a (i) of infinity. The reader is referred to [15, 19 for a thorough discussion of the 
details. 



4.1 The construction of C a 

Starting on the initial hypersurface <S, the construction introduced in [15] makes use of a blow-up 
of the point i £ S to the 2-sphere S 2 . This blow-up requires the introduction of a particular 
bundle of spin- frames over B a . Consider the (conformally rescaled) spacetime {M,g^) obtained 
as the development of the time symmetric initial data set (S, hij). Let SL(S) be the set of spin 
dyads S = {5a}a=o,i on S which are normalised with respect to the alternating spinor 6ab in 
such a way that eoi = 1. The set SL(S) has a natural bundle structure where S is the base space, 
and its structure group is given by SL(2, C), acting on SL(S) by S i— > <5 ■ t — {S A t A B }B=o,i- Let 
t = \/2eo, where eo is the future g-unit normal of S and 

t~aa' = g{r, S A ~5 A >) = e A °e A ,° + e A e A } 

is its spinorial counterpart — that is, r = r a e a = o- a AA ,T AA e a where c AA , denote the Inf eld-van 
der Waerden symbols and {e a }, a = 0, . . . , 3 is an orthonormal frame. The spinor taa> enables 
the introduction of space-spinors — sometimes also called SU(2) spinors, see (2j[7l[26]. It defines a 
sub-bundle SU(S) of SL(S) with structure group SU(2, C) and projection tt onto S. The spinor 
t aa allows to introduce spatial van der Waerden symbols via 

„AB (A B)A' a ^ A' a -i o q 

a a = a a A' T ' a AB= T ( B A)A'i & - M, <J. 

The latter satisfy 

h a b — CT aAB 0- AB , —5abO- AB 0- b CD — —£ A (c^D)B = h A BCD, 

with h a f, = h(e a , e^) = — 6 a b- The bundle SU(S) can be endowed with asu(2, C)-valued connection 
form uj a b compatible with the metric hij and 1-form <r AB , the solder form of SU(S). The solder 
form satisfies by construction 

h = h i jdx i ® dx j = h ABC D<y AB ® cr CD . (7) 



Given a spinorial dyad S £ SU (S) one can define an associated vector frame via e a = e a (S) = 
a AB S a t b b 5 B i, a = 1, 2, 3. We shall restrict our attention to dyads related to frames {ej}j=a,— ,3 
on B a such that e^ is tangent to the h- geodesies starting at i. Let H denote the horizontal vector 
field on SU(S) projecting to the radial vector 63. 

The fibre tt - 1 (z) c SU(S) (the fibre "over" i) can be parametrised by choosing a fixed dyad 
S* and then letting the group SU(2, C) act on it. Let (—a, a) 9 p 1— > 5(p 1 t A B ) S SU{S) be the 
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integral curve to the vector H satisfying 5(0, t A B ) = 5(t A B ) e tt With this notation one 

defines the set 

C a = {5(p,t A B ) g SU(B a ) | \p\ < a, t A B e SU(2,C)}, 

which is a smooth submanifold of SU(S) diffeomorphic to (—a, a) x SU(2,C). The vector field 
H is such that its integral curves through the fibre 7r _1 (i) project onto the geodesies through i. 
From here it follows that the projection map n of the bundle SU(S) maps C a into B a - 

Let, in the sequel, 1° = 7r _1 (i) = {p = 0} denote the fibre over i. It can be seen that 
1° w SU(2, C). On the other hand, for p e B a \ {i} it turns out that 7r _1 (p) consists of an orbit 
of U(l) for which p = \x(p)\, and another for which p = — \x(p)\, where x % (p) denote normal 
coordinates of the point p. In order to understand better the structure of the manifold C a it is 
useful to quotient out the effect of U(l). It turns out that X°/U(l) « § 2 . Hence, one has an 
extension of the physical manifold S by blowing up the point i to § 2 . 

The manifold C a inherits a number of structures from SU(S). In particular, the solder and 
connection forms can be pulled back to smooth 1-forms on C a satisfying the structure equations 
which relate them to the curvature form determined by the curvature spinor rABCDEF- In the 
conformally flat setting one has that 

rABCDEF = on B a . 



4.2 Calculus on C a 

In the sequel t A B G SU(2,C) and p <G R will be used as coordinates on C a . Consequently, one 
has that H = d p . Vector fields relative to the SU(2, C)-dependent part of the coordinates are 
obtained by looking at the basis of the (3-dimensional) Lie algebra su(2, C) given by 

0i\ 1/0 -1\ 1 / i 



Ul =2\i 0)' U2= 2\l ) ' U3 ^ 2 V -i 

In particular, the vector m 3 is the generator of U(l). Denote by Z i: i = 1,2,3 the Killing vectors 
generated on SU(S) by Ui and the action of SU(2, C). The vectors Zi are tangent to 1°. On J° 
one sets 

X+ = -(Z 2 + iZ 1 ), X_ = -(Z 2 -iZ 1 ), X = -2iZ 3 , (8) 

and extends these vector fields to the rest of C a by demanding them to commute with H = d p . 
It is noted that 

[X,X+]=2X+, [X,X_] = -2X_, [X + ,X_] = -X. 

The vector fields are complex conjugates of each other in the sense that for a given real-valued 
function W, X_W = X + W. More importantly, it can be seen that for p € B a \{i} the projections 
of the fields d p , X± span the tangent space at p. 

A frame cab — C(ab) dual to the solder forms <j cd is defined so that it does not pick 
components along the fibres — i.e. along the direction of X. These requirements imply 

(a AB ,c CD )=h AB CD , c CD = c 1 CD d p + c+ D X + +c CD X_, (9) 

where (•,•) denotes the action of a 1-form on a vector. In the conformally flat setting, from the 
properties of the solder form o AB one finds that 

1 1 

C AB = X AB , c\ B = -ZAB, C AB = -VAB, (10) 

with constant spinors xab, Dab and Zab given by 

_ Ai 1 lii loo 

X A B = V2e (A € B) , UAB = ^^ e A e B ) Z AB = ~^ € A e B ■ 
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The connection coefficients, jabcd, are defined by contracting the connection form with the 
frame cab- In the conformally flat case, one has 

1ABCD = 7T - {^AC X BD + (BD^AC)- (H) 

2p 

Let / be a smooth function on C a 

Dab! = cab (f)- 

Similarly, let [Iab be a spinor field on C a . Then the covariant derivative of /jab is given by 

Dab^cd — cab(^cd) — 1ab E cI i ed — 1ab E d^ce- 
Analogous formulae hold for higher valence spinors. 

4.3 Normal expansions at % 

In [15] a certain type of expansions of analytic fields near i has been discussed. Suppose 
£a b —A b denotes the components of an analytic even rank spinorial field with respect to a 
fiduciary spin frame 5^- One can introduce the expansion 

Cmb^bM) = E -^\ Pn ° pDp ■ -n ClDl D Cp D p ■ ■■Dc 1 d 1 Ca 1 b 1 ...a 1 bM ( 12 ) 

with n AB = n AB (q), q 6 B a , the spinorial representation of the vector n l di — (x 1 /\x\)di, riin 1 = 
-1. 

4.4 An orthonormal basis for functions on SU(2,C) 



The lift of the expansion (|12[) from B a to C a introduces in a natural way a class of functions 
associated with unitary representations of ST/ (2, C). Namely, given t A B e SU(2,C), define 

/ \ 1/2 / \ 1/2 



with j,k = 0, . . . , m and m = 1, 2, 3, . . .. The subindex expression (Ai---A m ) k means that the 
indices are symmetrised and then k of them are set equal to 1, while the remaining ones are set 
to 0. Details about the properties of these functions can be found in [Til HI]- The functions 
\fm + VT m 3 k form a complete orthonormal set in the Hilbert space L 2 (fi, SU(2, C)), where \i 
denotes the normalised Haar measure on SU{2, C). In particular, any analytic complex-valued 
function / on SU(2,C) admits an expansion 

oo m m 

f(t A B) = ^ ^ fm,k,j T m fc 7 ( tA B ) ) 
m=0 j=0 fe=0 

with complex coefficients f m ,k,j- Under complex conjugation the functions transform as 



Tj k = (-iy +k T m m -\ n _ k . (13) 
The action of the differential operators (JH) on the functions T m is given by 



XT m k j = (m-2j)T m k j , 

X+T m k 3 = y/jim-j + VT^, X_T TO * = -v/(j + l)(m-j)T m fe : 
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A function / is said to have spin weight s if 

Xf = 2s f. 

Such a function has an expansion of the form 

oo m 

/ = ^""^ fm,kT m m /2-s- 

m>|2s fe=0 

Finally it is noted that products T t Jl k x T i2 ■"^ can be linearised — that is, written as a linear 
combination of other functions T. J , , for suitable i, j, k using the formula 

P 

T ii J1 fe 1 x T i 2 h k2 =^2 D ( i i'h, k i; i 2,j2,k 2 ;ii + i 2 - 2p, ji + j 2 -p, fci + fc 2 -p) 

with p = min{ii, i 2 ,ji + j 2 , &i + fe} and 



D(ix,ji,ki;i 2 ,h, l,m,n) = C(ii,ji;i 2 ,j 2 ; l,m)C(ii,ki;i 2 , k 2 ;l,n), 
and C(i,j; k, I; m, n) the Clebsch-Gordan coefficients of SU(2, C). 



4.5 Normal expansions at X° 

In the sequel, it will be necessary to lift analytic fields defined on B a to C a . As in section |4~3I 
consider normal coordinates x % on B a centred on i and which are based on the orthonormal frame 
c *a = a a c *AB ~ a a B ^A^B- ^ n terms of p and t A B on C a and the normal coordinates x l , the 
projection ir has the local expression 

it : (p,t A B ) - - v^^o^i- 

This expression can be used to carry out lifts to C a . In particular, the lift of \x\ is p. Applying the 
procedure described in [15j to the expansion (|12p one obtains the expansion of the spinor- valued 
function £a 1 b 1 ---AiBi on C a . Denote by £j = £,(A 1 B 1 ---AiBi) j i < j < I its essential components. 
The function £j has spin weight s = I — j and a unique expansion of the form 

00 

^ = Ew, (15) 

with 

p+l 2 q 

^ \P = X! £j.P;2g,fc ^2(7 2o-j+7' ( 16 ) 

g— max{|Z — p} k—0 

and complex coefficients £j p - 2 q,k- More generally, we shall consider symmetric spinorial fields 
£,Ai—A r on Ca with independent components £j — £(A 1 ---A 2r ) ■■> < j < 2r, and spin-weight 
s = r — j which do not descend to analytic spinor fields on B a - In this case one has that 

00 q(p) 2q 

p=0 q=|i — j| fe=0 

where one has a priori that < |r — j| < q(p). An expansion of the latter form will be said to be 
of type q(p). 
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4.5.1 On the expansions of W on C a 

From the previous discussion it follows that the function W on B a admits a lift to C a . This lift 
will be again denoted by W . Its normal expansion at 2° is given by 

oo 2p 
p=2 k=0 y ' 

with w p -2 Pt k S C given by 

w p , 2p , k = (V2Y Q ^ 1/2 D (BpCp ■ ■ ■ D BlCl)h W{i), 

=(^) p (?) 1/2 ( 2 ;)" 1 % 1 c I --u )fe ^---v 

Thus, Wi 1 ...i p = if and only if w p . Pt k = 0, k = 0, . . . , 2p. The function IU is the lift to C a of a 
real function. Hence, it satisfies W = W . Using the property (fT5|) it follows that the coefficients 
w P ;2p,k satisfy the reality condition 

Wp- 2p ,k = (-l) p+fc w p;2 p,fc, ft = 0, . . . ,2p. 

In particular, the coefficients w p -2p t p are real. If the function W is axially symmetric, then its 
normal expansions on C a take the simpler form 

OO -. 

TTL V — ^ t r> 71 

5 The spacetime Friedrich gauge 

The formulation of the initial value problem near spatial infinity presented in [15] employs gauge 
conditions based on timelike conformal geodesies. The conformal geodesies are curves which 
are autoparallel with respect to a Weyl connection — i.e. a torsion-free connection which is not 
necessarily the Levi-Civita connection of a metric. An analysis of Weyl connections in the context 
of the conformal field equations has been given in [14] . In terms of this gauge based on conformal 
geodesies — which shall be called the Friedrich gauge or F-gauge for short — the conformal factor 
of the spacetime can be determined explicitly in terms of the initial data for the Einstein vacuum 
equations. Hence, provided that the congruence of conformal geodesies and the fields describing 
the gravitational field extend in a regular manner to null infinity, one has complete control over 
the location of null infinity. This can be ensured by making B a suitably small. In addition, the 
F-gauge renders a particularly simple representation of the propagation equations. Using this 
framework, the singular initial value problem at spatial infinity can be reformulated into another 
problem where null infinity is represented by an explicitly known hypersurface and where the data 
are regular at spacelike infinity. The construction of the bundle manifold C a and the blowing up 
of the point i £ B a to the set 1° C C a , briefly described in section [4] are the first steps in the 
construction of this regular setting. The next step is to introduce a rescaling of the frame bundle 
so that fields that are singular at 1° become regular. This construction is briefly summarised in 
this section. 



5.1 The manifold M a , K 

Following the discussion of 15] assume that given the development of data prescribed on B a , 
the timelike spinor t aa introduced in section 3] is tangent to a congruence of timelike conformal 
geodesies which are orthogonal to B a . The canonical factor rendered by this congruence of 
conformal geodesies is given in terms of an affine parameter r of the conformal geodesies by 

e = n^n (i - with u = - 2n =, (17) 



n 



where fl — i9~ 2 and i3 solves the Yamabe equation ([3]) — see [14j [X5l [E]. The function k > 
expresses the remaining conformal freedom in the construction. It will be taken to be of the form 
k = k'p, with k' analytic, ~ 1. Associated to the conformal factor there is a 1-form d p 
from which the Weyl connection can be obtained. In spinorial terms one has that 

Aaa 1 — -^J=TAA'd T Q - T B A ,d,AB, 

where cLab is calculated in the case of conformally flat data via 

'x A b - p 2 D AB W 



d A B = 2p 



(1 + P Wf 



The function k in the conformal factor 0, induces a scaling 5a i_ * k 1 ^ 2 $a of the spin frame. 
Accordingly, one considers the bundle manifold C a , K = K 1//2 C a of scaled spinor frames. Using C a K 
one defines the set 

M a , K = Ur^)\ q eC a , K ,-^-<r< Uj{q) 



n(q) n(q) 

which, assuming that the congruence of null geodesies and the relevant fields extend adequately, 
can be identified with the development of B a up to null infinity — that is, the region of spacetime 
near null and spatial infinity. In addition, one defines the sets: 

!={(T,q)eM a , K \p(q) = 0, M<1}, 
Z ± = {(r,q) e M a , K | p(q) = 0, r = ±l}, 

^ ± = f (r, g) € M a , K | p(q) > 0, r = ±^4 

which will be referred to as, respectively, the cylinder at spatial infinity, the critical sets and 
future and past null infinity. In order to coordinatise the hypersurfaces of constant parameter r, 
one extends the coordinates (p, t A B ) off C a>K by requiring them to be constant along the conformal 
geodesies — i.e. one has a system of conformal Gaussian coordinates. 

Remark. For the purposes of the analysis carried out in this article it turns out that the most 
convenient choice of the function k in the conformal factor of equation (|17[) is 

k = p. 

This leads to considerable simplifications in all the relevant expressions. From this point onwards, 
this choice will always be assumed. 



5.2 The conformal propagation equations 

On the manifold A4 a , K it is possible to introduce a calculus based on the derivatives d T and d p 
and on the operators X+, X~ and X. The operators d p , X+, X^ and X originally defined on 
C a can be suitably extended to the rest of the manifold by requiring them to commute with the 
vector field d T . In order to derive the propagation equations, a frame caa' and the associated spin 
connection coefficients Taa'bc of the Weyl connection V will be used. The gravitational field is, 
in addition, described by the spinorial counterparts of the Ricci tensor of the Weyl connection, 
®AA'BB<, and of the rescaled Weyl tensor, <Pabcd — see [HI [15j [19] . Some further notation will 
be required. Let 

(f>ABCD = 4>oe°ABCD + <l>l e ABCD + fa^ABCD + fo e ABCD + <t>± e ABCDi 

where 

4>i = 4>(ABCD)i, e ABCD = e (A^ Ee B F ' e C Ge D) H '' 1 « = 0, . . . , 4, 

where expressions like (ABCD)i mean that after symmetrisation i indices are set to 1. A space 
spinor Qabcd = Qab(cd) is introduced such that 

Qaa'cc — Qabcdt B ait d c < 
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The space spinor Qabcd shall be further decomposed as 



QABCD — Q(AB)CD + -Z e AB@G CD- 



From the spin coefficients V aa' bc one defines 



Fabcd — t b b TaB'CD, 
which in turn, will be decomposed as 

Fabcd — (£abcd — X(ab)cd) — ^abIcd- 
They possess the following symmetries 

FaBCD = ^AB(CD), XABCD = XAB(CD), £,ABCD = £,(AB)(CD)- 

The term £,abcd is related to the intrinsic connection of the leaves of the foliation defined by 
taa 1 ; the spinor X(AB)cd corresponds to the second fundamental form of the leaves; and /ab is 
the acceleration of the foliation. 

Using the F-gauge it can be shown that the extended conformal field equations given in 
15J imply the following evolution equations for the unknowns c AB [p, — 0, 1,±), £abcd. Jab, 

X{AB)CD, Q(AB)CD, ©g^H 1 

EF „0 



drCAB = -X(AB) EFc EFi ( 19b ) 



9rC°AB = -X(AB) C EF ~ fABi (19a) 
J 

C (AB) 

d T £,ABCD = -X(AB) EF ^ EF CD + -^(^ACX(BD)EF + ^BDX(AC)EF)f EF 

-^X(AB)(C E ^ fD)E - g (£AC<d F F BD + e BD<d F F Ac) ~ i®V>ABCD, (19c) 

drfAB = -X {A B) EF fEF + ^=Q F F A B, (!9d) 

drX(AB)CD = -X(AB) EF ^ EF CD ~ &(CD)AB + QllABCD, (19e) 
d T <d(AB)CD = -X( C D) EF ®(AB)EF ~ d T QrjABC D + iV2d E \ A HB)CDE, (19f) 
dr& G G AB = -X ( AB) EF ®G G EF + V2d EF VA BEF, (19g) 

where 

VABCD = -^{4>ABCD + T A A Tg B T c ° T D D <j) A > B'C> D>), 
HABCD = —z{<t>ABCD - T A A T B B T C C T D ° ^ A 'B'C'D')> 

denote, respectively the electric and magnetic part of of <\>abcd- Thus, the equations (|19ap - |T9g| 
are essentially ordinary differential equations for the components of c AB , Tabcd, Qabcd- The 
redundancies in the latter equations, which are due to the symmetries of the spinors, can be 
eliminated by noting that valencc-2 spinors can be written as a linear combination of the spinors 

AT 1 _ 1 i i _loo 

X A B = V2e (A € B) , y A B = ~^ € A e B i Z AB = ~J=t . A ^B ) 

while the valence-4 spinors can be written in terms of the spinors £abcd (* = 0, . . . , 4) , Habcd, 
xaccbd + xbd£ac, Vaczbd + Vbd^ac and zac^bd + zbdcac- The evolution equations for the 
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spinor 4>abcd are derived from the Bianchi equations. One has the following Bianchi propagation 
equations: 

(y/2 - 2c o , 1 )d T o + 2c° ^^ - 2dg 1 d a <h + 2cg,0 o 0i 

= (2rooii - 8rioio)<^o + (4r 00 oi + 8riooo)</>i - 6r ooo^2, (20a) 
V29 r 0i - c^drcpo + CQ O a T 2 - cfidafa + c% d a <j)2 

= -(4Tmo + /n)0o + (2r on + 4rnoo - 2/oi)^i + 3,/oo02 - 2r OO oo03, (20b) 

= -rim^o - 2(r 110 i + /n)0i + 3(r 00 n + rnoo)02 

-2(r 001 ~ /oo)03 ~ 1^00004, (20c) 
V29r03 - C?i9r02 + Co O d T 04 ~ cji9 a 02 + CoO<9«04 

= -2rmi0i - 3/n02 + (2rnoo + 4r 00 ii + 2/oi)03 - (4r 00 oi - /oo)04, (20d) 
(V2 + 2c o i 1 )gV0 4 - 2c? 1 a r & + 2c o " i a Q 4 - 2c? 1 0«& 

= -6r im 2 + (4r mo + 8r om )03 + (2r U00 - 8r O i O i)04, (20e) 

with a = 1,± and d\ = d p , d± = X±. In addition to the latter propagation equations, we shall 
also make use of a set of three equations, also implied by the Bianchi identities, which we refer 
to as the Bianchi constraint equations: 

Cn<9 T 0o - 2(% l d T <j>i + cj] o <9 T 2 + c^d^o - 2c o \<9 Q! 0i + Co O <9 Q 2 
= — (2r (01)11 — 4r m o)0o + (2r on — 4T( i)oi - 4rn OO )0i 

+6r( O i)oo02 - 2r OO oo03, (21a) 
Cn<9 T 0i - 2c^ 1 S T 2 + Cq O <9 t 3 + c" 1 5 a 0i - 2c o \<9 Q! 02 + Co O <9 a 3 

= rmx0o — (4r (01)11 — 2rnoi)0i + 3(roon — rnoo)02 

-(2r 00 oi - 4r( O i)oo)03 - r OO oo04, (21b) 

c?i9 T 02 - 2c° l( 9 T 03 + c° o <9 T 4 + Cii<9 Q 02 - 2co 1 <9 Q 03 + Co O <9 Q 4 

= 2r m i0i - 6r( O1 ) U 2 + (4r on + 4r (0 i)oi - 2rn O o)03 

-(4r ooi - 2r(oi)oo)04- (21c) 

The propagation equations are supplemented by initial data on C a K constructed from the 
conformal factor il — d~ 2 and the flat connection coefficients Jabcd — cfr. equation (fTTj) — by 
using the conformal constraint equations — see e.g. |14j for more details. One finds that for 
conformally flat data and the gauge choice n = p: 



®abcd — ^-D(ab-Ocii)^i (22a) 
P 3 

4>abcd = -q2-D(abDcd)Q, (22b) 

c% = 0, (22c) 

C AB = PXAB, (22d) 

cab = z ab, (22e) 

cab = Vab, (22f) 

Ubcd = 0, (22g) 

X(AB)CD = 0, (22h) 

fAB=X A B, (22i) 



on C a ,K- 
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5.3 Structural properties of the evolution equations 

We discuss now some general structural properties of the equations ()19a[) - p^g| ) , (|20b[) - (|20ep and 
()21a|l - (|21c|) which will be used systematically in the sequel. Introduce the notation 

v = {c°ABi c AB^ABCD-, ®ABCd) , 4> = (00, 01,02,03,04) ■ 

The unknown vector v has 45 independent complex components, while has 5 independent 
complex components. In terms of v and 0, the propagation equations ()19ap -( [T9gl ) can be written 
as: 

d T v = Kv + Q(v, v) + L<f>, (23) 

where K and Q denote, respectively, a linear and a quadratic constant matrix-valued function 
with constant entries and L is a linear matrix-valued function with coefficients depending on the 
coordinates and such that L\ p= q = 0. Similarly, the systems (|20bjl - (120cp and ()21all - (121c|l can be 
written as 

x/2£cV0 + A AB c AB d^ = B{T ABCD )<t>, (24a) 
F AB CA B d^ = H(T ABCD ) : (24b) 

where E denotes the 5x5 unit matrix and A AB c AB , /i = 0, . . . , 3, are 5x5 matrices depend- 
ing on the coordinates, while B(T abcd) denotes a constant matrix- valued linear function of the 
connection coefficients Tabcd- On the other hand, F AB c AB denote 3x5 matrices with coordi- 
nate dependent entries and H (T abc d) is another constant matrix- valued linear function of the 
connection coefficients T AB cd- 

Consider now the system (l2"3" |) - (|24ap with data given on C a . R . Given a neighbourhood W of 
C a , K in M a .K on which a unique smooth solution of the Cauchy problem is given. From the point 
of view of the propagation equations, the subset Wnl is a regular hypersurface. Introduce the 
notation = v\wni, 4>^ = 0|wnx- Due to the property £| p= o = 0, equations (|23|) decouple 
from equations (|24a|) and can be integrated onWfll Using the observation that the restriction 
of the initial data to 1° coincides with Minkowski data one has that on Wnl: 

Q (0) =0, xfl B)CD =^ f { AB=XAB, i { A ] BCD = ^ (^) 

(c° AB )^=-rx AB , (4 B ) (0) =0, {c- A J 0) =y A B, {c + AB r^z A B. (25b) 
It follows that A 1 = A AB c\ B is such that 

A |wnx = 0, 

so that the system (|24al) also implies an interior system on Wfll whose solution is 

4> { ABCD = ABCD- (26) 

The solutions (|25a)) - (I25b|) and P5|) extend analytically to the whole of X and in particular to the 
critical sets T ± . The set X is a total characteristic of the system (|23|) - (|24a.|) in the sense that the 
whole system reduces to an interior system on X. Moreover, the constraint equations (|24b[) also 
reduce to an interior system on X. Another crucial structural property is that 

A = V2E + A AB c° AB = \/2diag(l + t, 1, 1,1,1- r) on X, 

so that the matrix ^4° which is positive definite degenerates at X^ . Understanding the effects of 
this degeneracy is the main motivation behind the analysis in this article. 

The previous discussion can be generalised by repeated application of the differential operator 
dp to the equations (f2"3")l . (|24a[) and (|24b[) to obtain interior systems for the quantities = 
dp P ^v\j and 0^ = 9p P V|i which will be called the order p transport equations. Their behaviour 
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on the whole of X will be studied in the sequel. The transport equations then take the following 
form for p > 1 : 

d T v^ = Kv^ + Q(w (0) ,f (p) ) + QO (p) y o) ) 

+ ( P ) (Q( v(j) ,v (p ~ j) ) + L^^P-^j +L^(j) (0) , (27a) 
(V2E + A AB (c%)^)d T ^ + A AB {c%)^d^ = B{Yf BCD )^ 

+ E ft) { B ^abcd)4> { ^ ] - A AB (c%)^d^ ] ) , (27b) 
F^(c%)(°)a r ^) +F^(c^)(°)5^W = H{Tf BCD )^) 

+ E ft) (^'L^^ - F AB {c p AB )^d^) . (27c) 

In the previous equations the values of i/ * 1 and cj>^ given in (125a[) - (|25bp and (|2"6")l are assumed. 
Note that the non-homogeneous terms in the equations (|27ap - (|27c[) depend on v^ p \ <p i - p ^ for 
< p' < p. Thus, if their values are known, then (|27a|) - (I27b|) constitutes an interior system of 
linear equations for and <fi^ ■ The principal part of these equations is universal, in the sense 
that it is independent of the value of p. If the initial data on C a , K for the system (f2Ul)- (|24aj) is 
analytic — as it is the case in the present analysis — then suitable initial data for the interior 
system (|27a[) - (|27bj) can be obtained by repeated ^-differentiation and evaluation on J°. 

The interior system (|27aj) - (j27b[) is decoupled in the following sense: if v^ p \ 4>^ p > are known 
for < p' < p one can solve first (|27a[) as it contains at most quantities of order With 
the knowledge of v^' at hand one can then solve (|27b[) to obtain </>( p ). 

The language of jets is natural in the present context. For p = 0, 1,2, . . . and any sufficiently 
smooth (possibly vector valued) function / defined on M a , K , the sets of functions {f^°\ ■ ■ ■ , 
y(p)} on x will be denoted by [/] and referred to as the jet order p of f on X — and similarly 
with X replaced by X°. If u = (v, </>) is a solution to the equations (|27a|) . (|27b() and (|27cp . we refer 
to 4 p) [u] as to the s-jet of u of order p and to the data [u] as to the d-jet of u of order p. An 

s-jet J^jf 1 [u] of order p will be called regular on X = X U X + U X~ if the corresponding functions 
extend smoothly to the critical sets X ± . 



6 The Schwarzschild spacetime in the F-gauge 

Due to its relevance for our purposes, a brief discussion of the Schwarzschild spacetime in the 
F-gauge is provided. The material is adapted from the one given in [15] . The Schwarzschild line 
element with mass m in isotropic coordinates is given by 

g= '-j-) dt 2 -(l + — (dr 2 + r 2 da 2 ) , 

\1 + m/2r J \ 2rJ K ' 

where da 2 = d6 2 + sin 2 9dip 2 is the standard line element of the unit sphere § 2 in polar coordi- 
nates. Writing the first fundamental form h, t j and the second fundamental form Xij on the initial 
hypersurface 5* = {t = 0} in terms of the coordinate p — l/r one obtains 

hij ^ hij , Xij 0, 

where 

h = hijdx'dat = - (dp 2 + p 2 da 2 ) , Q = (1 + ^ /2)2 . 
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Thus, we have initial data set for the Schwarzschild spacetime which is time symmetric and 
conformally flat. Accordingly, one has that 

m 

U=l, W = -, 

near p = 0. In what follows let k = p and assume that a is chosen small enough such that 
(1 + (m/2)p) 7^ for \p\ < a. The conformal factor 8 and the 1-form cLab associated to the 
F-gauge read in this case 

H . " .J \ . " A d AB = 2PXAB 



(l + pm/2) 2 V (l + pm/2)V' (l + pm/2) 3 ' 

Furthermore, the non-trivial initial data on C a , K is given by 



_ 6mp 2 G 
°(ab)cd- {1 + pm/2)2 e ABCD , O GC j> 

,2 



4>abcd = -6me ABCD 

This initial data set for the propagation equations (|19ap - flljg| ) and ()20a[) - (|20e[) is explicitly spher- 
ically symmetric — the functions involved are of spin-weight and contain only the harmonic 
T ° . Accordingly, one puts forward a spherically symmetric Ansatz for their development. More 
precisely, one writes 

C AB = c x x AB, CAB — c x x AB, C AB — C y IJAB, c\b = c t z AB, 
fAB = fxXAB, £,ABCD = ^xi^ACXBD + <^BD<^Ac)i 
X(AB)CD = X2 e ABCD + Xhh-ABCD, 

®(AB)CD — ^>2^ABCD + ©ft^ABCU, ©G CD = @x%AB i 
4>ABCD = 4>2£ 2 a BC D> 



where 



c\, c+, c , f x , £ x , £ 2 , Xh, ©2, ®h, 



are functions depending on (p, r) only. The previous Ansatz, together with the propagation 
equations p9ap - |T5g| and (|20a|) - (|20cp imply an initial value problem of the type 

d T u — F(u,t, p;m), u(0, p; m) = uo(p; m), (28) 

with analytic functions F and uq for the unknowns 

U = { c x, cl, C~ , C+ , fx, Xx, £.2, £,h, 92,Qh, 6x, fo) ■ 

The solution with m = corresponds to a portion of the conformal Minkowski spacetime in which 
the only non-vanishing components of the solution are given by 

c x = — r, c x = p 1 = Cy = 1, f x = 1. 

Since in this case the solution exists for all r, pel, it follows from standard results of ordinary 
differential equations that for a given m there is a sufficiently small po such that there exists 
an analytic solution to the system (|2"5)) which extends beyond for \p\ < \pq. Hence, one can 
recover the portion of the Schwarzschild spacetime which lies near null and spatial infinity if a is 
taken to be small enough. 

It follows from the above discussion that the coefficients that are obtained from solving the 
transport propagation equations on the cylinder at spatial infinity correspond to the terms in the 
Taylor-like expansions 

u = E 
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of the solutions of the initial value problem l|28j) . In particular, the logarithmic singularities that 
have been observed to appear in the critical sets of the development of more general classes 
of data — see e.g. [3D] — do not arise in the case of the Schwarzschild spacetimq^J. One has the 
following result. 



Proposition 1. The solutions of the transport equations \27a\ ) and {27b) for time symmetric 
Schwarzschild initial data extend analytically through T^- for all orders p. Moreover, the solutions 
to the transport equations are polynomial in t. 

For the purposes of the present article it turns out that it will be necessary to know the 
expansions explicitly up to order p — 4 (inclusive). These straightforward, but nevertheless 
lengthy computations have been performed with the aid of a computer algebra system Maple V. 

7 Further properties of the Bianchi transport equations 

The transport propagation equations for the Bianchi subsystem of interior equations (|27bj) read 
explicitly 

(29a) 
(29b) 

d T ^> + -X + ^> + -X^T = RT, (29c) 
dA P) + + = B!f\ (29d) 

(29e) 

with Rj — Rj(u^°\ . . . , u( p_1 )), j = 0, ... ,4. On the other hand, the Bianchi transport constraint 
equations (|27cp are given by 

(30a) 
(30b) 
(30c) 

with Sj = Sj(u(°\ u( p -V), j = 1, ... ,3. 
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In order to extract detailed information from the above equations one makes use of an explicit 

i k ' 



decomposition of the various functions in terms of the spherical harmonics T- J ,. We recall the 



following lemma which was proved in [T5] . 

Lemma 1. The following rules for expansion types hold: 

(i) The functions (c\ B — pxabY p \ v^ p \ <t^ p \ p — 1,2, .. . on T are of expansion type p — 2, 
p — 1, p respectively. 

(ii) The functions R!f\ i = 0, ...,4 and Sj, j — 1,2,3 are of expansion type p—1 for 
I' 1-2 



(Hi) If for a given integer p > 1 the data for <fi( p ' on C a ^ K are of type p — 1, then 4>^ p ' on T is of 
type p—1. 



iMore generally, this has been shown to be the case for all asymptotically flat vacuum static spacetimes — see 

EH. 
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7.1 Decomposition in terms of spherical harmonics 

Given the vector = (iq , . . . ,1$) — respectively «W, <f>^ — and non-negative integers q 
and k = 0, . . . , 2q one defines the sector & q ^[u^\ as the collection of coefficients 



Ui; 2q ,k = (2g + 1) / u ( f ] T 2 fe d/i, 
Jsu(2) 



5(7(2) 

where s is the spin- weight of uf\ and d/i is the Haar measure of SU(2). Furthermore, one defines 

2q 

& q [vF>] = (J 6,,* [«<">]. 

A sector will be said to be vanishing if & q [u^] — {0}. 

The Weyl spinor of time symmetric, conformally flat initial data is of expansion type p — 1 on 
C a , K - Accordingly, one writes 



P 2q -1 



g=|2-i| fe=0 



X/ X/d) a j,p;2q,kT 2q q _ 2+ j> 



with complex (r-dependent) coefficients dj tP -2 q ,k- The normalisation factor ( .) has been added 
for convenience. The substitution of the latter expression into equations (|29ap - (|29e[) and (I30a,p - 
(|29ap renders equations for the various coefficients o?,p ; 2g,fc- In the cases p > 0, q = 0, one finds 
the equations 

and 

ra 2,p;0,0 — P a 2,p;0,0 = 6S2,p;0,0- 

If p > 1, <? = 1, fe = 0, 1, 2 one finds 

a 'l,p;2,fc + 2^202, p;2,fc + Ol,p;2,k = 4i?l,p ; 2,fc, 



3 3 

a 2,p;2,fc + ^/^2 a 3,p;2,fc — ^ P2«'l,p;2,k = 6i?2,p;2,fc, 
1 



°3 : p;2,fc — g/32a2,p;2,fc — a 3,p:2,k — 4i?3, p; 2,fc, 

and 

3 3 

Ta 2,p;2,fc + ^/?2fl3,p;2 : fc + ^/^O^p^fe ~ Pa2,p;2,fc = 6S2,p;2,fc- 

More crucially, one obtains for 2 < p, 2 < fe = 0, . . . , 2q the equations 



and 



(1 + T)a' Q p . 2q k + ^PiahpVq* - - 


2)ao,p;2tj,fc = 


Ro,p;2q,ki 


(31a) 


a l,p:2q,k + 


g/^2fl2, p;2<j,fc — 2/3iao,p ; 2 (? ,fe 


+ «l,p;2g,fc = 


4Rl,p;2q,ki 


(31b) 


a 2,p;2q,k + 


-^02O-3.,p;2q,k ~ ^/?2 a l,p;2g,fc 


= 6i?2,p;2q,fe 




(31c) 


a 3,p;2q,k + 


2/3ia4,p;2(j,fe — g/32 a 2,p;2q,fc 


~~ a 3,p;2q,k — 


4^?3,p:2g,fe: 


(31d) 


(1 - r)ai, 


»;2<j,fc ~ -^Pl a 3,p;2q,k + (P ~ 


2)a4,p;2ij,/c = 


-^4,p;2ij,fc) 


(31e) 


Ta 'l,p;2q,k + 


-^020-2, p:2q.k + 2/3iao,p;2q,fc 


— PO>l,p;2q,k - 


= 45i ]P; 2g,fc, 


(32a) 


Ta 2,p:2q,k + 


-^P2<l3,p\2q,k + 4 P2«l,p;2g,fc 


— P(l2,p:2q,k 


= 652,p;2g,/s; 


(32b) 


Ta 3,p:2q,k + 


2/3ia4 jP; 2 9l fc + g/?2a2,p;2q,fe 


— Pd3,p;2q,k ' 


= 4S'3 ) p;2 g ,fc, 


(32c) 
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where 

ft = V(g-i)(g + 2), /3 2 = Vg(gTl), 

and Ri tP -2q,ki j = 0, . . . ,4 and Sj tP -2 q ,k, i = 1,2,3 are such that 

g=|2-j| fc=0 g=|2-j| k=0 

The functions i?| p) and sj p) contain products of 0( p ') and tA") for < p' < p - 1, < p" < p - 1 
so that in order to obtain the representation (I33|) one has to linearise products of the form 
T. j \ xT. h . using the formula (HI). 

For latter reference the following result is noted. 

Lemma 2. // the s-jets j!f ^ [v] and j!f ^ [(/)] have polynomial dependence in r for some p > 1, 
then J^[v] has also polynomial dependence in r. 

The proof of this lemma follows directly from the structure of the transport equation (|27a[) . 
7.2 Discrete symmetries of the development 

It is well known that if a spacetime {M.,gfiv) is the development of time symmetric data, then 
the spacetime has a discrete time reflexion symmetry: that is one has g^ftyX ) = g^ v {— t, x a ), 
where t is a time function such that t = yields the slice of time symmetry. However, in order to 
discuss the effect of this discrete symmetry on spinorial objects one has to be more careful as the 
transformation t i— > —t changes the handedness of the canonical orthonormal tetrad associated 
to a spin dyad {Sa}a=q.i- 

In the case of the manifold M a ,K, the discrete transformation t t— > —t induces on the dyad 
{5a}a=o,i the transitions 

S i-> Si, Sx i-> <5 , 
so that, for example, taa 1 l— * taa 1 ■ On the other hand, one has that 

£ab !-* — cab, 

XAB !-» XAB, VAB ^ Z AB , Z A B >-> "VAB- 

Similarly, one has that 

e ABCD l— * e ABCD' i = 0, . . . , 4, 

(eiciBD + (bdxac) >-> -(f/tcSBD + (bdxac), 

(CACVBD + tBDVAc) (eAC^BD + Csn^Ac), (eACZBD + fBDZAc) !-> {tACVBD + ^bdVac), 
h-ABCD | — > hABCD- 

From the latter, one deduces the following transition rules on the components of the Weyl spinor 
0abcd: 

4 (t) i-> o (-t), 03 (r) h-> 0i(-t), 2 (r) 2 (-r), 0i(t) h-> 3 (-t), o (t) h-> 4 (-t). 

Furthermore, in order to correct the change of handedness produced by the transformation r i— > 
— r, one has the following correspondence rules for the operators X±: 

X+^-X^, X— i — > —X+. 

In addition, it is noted that 

9 T 4 (r) i ^ -St0o(-t), 9 r 3 (r) h-> -<9 t 0i(-t), 9 T 2 (r) i-> -9 r 2 (-r). 
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Combining the above rules together with equations (|29ap - (|29ej) and (|30ap - (|30cp one deduces the 
following transition rules for their right hand sides: 

Ro >-> -Ri, Ri i-> -R3, R2 ^ -R2, R3 ^ ~Ri, Ri >-> -Ro 
Si 1-* 5*3, S 2 i-> 52, S 3 i-> Si. 



For a spacetime with time reflexion symmetry, the aforediscussed transition rules allow to 
deduce the symmetries satisfied by the various terms appearing in the transport equations (|20a[) - 
pOe] and ([2la |) -([2Tc ]) . 



Lemma 3. For a spacetime arising from time symmetric initial data, the solutions to the Bianchi 
transport equations satisfy the parity conditions 

^a 3 ^ ;29: fc(r) = (-l) r ^ai iP;2( j,fc(-T), 

^r a 2,p;2g,fe(r) = (- l) r ^a2,p ;2? ,fc (~t) , 

and afeo 

a;i? 4 , P ;2 9 , fc (r) = (-iy +1 d r T R , P ;2 q A-r), 

d r T R 3 , p;2g ,k(T) = (-l) r+1 d r T R hP ;2 q A-r), 
d r T R 2 , M {r) = (-l) r+1 5;i? 2 , p;2g , fe (-r), 
d r T S M , k (r) = (-l) r d r T S hp , 2q . k (-r), 
5;5 2 , p; 2 9 ,fc(r) - {-l) r d r T S 2 , p , 2q . k {-r). 



7.3 A procedure to solve the Bianchi transport equations 

A first analysis of the structure of the solutions to equations (|31aj) - (l31ej) has been given in [T5] . 
In particular, in the aforementioned reference a procedure was given by means of which the 
constraint equations (|32a|l - (|32cp are used to eliminate the unknowns ai, p; 2g,fe, a2, P -,2 q ,k, a<3, P ;2 q ,k 
so that to find a solution to the transport equations (|31ap - (|31c[) it is only necessary to solve a 
reduced system involving ao jP]2 q,k and a^ p - t2q ^- The remaining coefficients are then obtained by 
means of purely algebraic manipulations. 

For the purpose of the present investigation it will turn out to be more convenient to consider 
an alternative procedure to find solutions of the equations (|31al) - (|31el) . Again, equations (I32a|) - 
(I32cp will be used to obtain a reduced system. But in this case, the system will involve a\^ 2q ^ 
and a3 iP; 2g,fc. Due to the formal similarity with the reduced systems obtained in [34], this type of 
reduced system will be called Maxwell-like. The reasons to prefer this approach over the one put 
forward in [15] will be explained towards the end of this section. 

In the following discussion, the values of the indices p, q and k are considered as fixed. In 
order to ease the formulae, obvious indices will be suppressed in the following. 

One can construct a MaxwelHike propagation system involving the coefficients 01, a 2 and 
03 by considering the sum of equations (|31bp and (|32a|) . equation (|31cp and the difference of 
equations (|31dp and (|32c|) . The resulting equations are given by 

(1 + r)oi + ^(3 2 a 2 - (p - l)oi = 2Ri + 45i , (34a) 

a 2 + ^/3 2 a 3 ~ \P2a1 = 6R2, (34b) 

(1 - r)a 3 - 2 -(i 2 a 2 + (p - l)o 3 - 2i? 3 - 45 3 , (34c) 

where the prime denotes differentiation with respect to r. Note that the above equations do not 
contain oq or 04. The associated Maxwell- like constraint equation is given simply by equation 
(j32b|) . Namely 

3 3 

™2 + 7/^203 + jfoai - pa 2 = 6S 2 . (35) 
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From equations ([35| and (|34b[) one obtains the algebraic relation 



3 3 

-/3 2 (1 - r)a 3 + -/3 2 (1 + r)ai -pa 2 = 6S 2 - 6ri? 2 , 



(36) 



which can be used, in turn, to eliminate a 2 from both (|34a[) and (|34cp so that 
(1 + r)oi + - p + ^g(<? + !)(! + r)j a x + -^g(a + !)(!- r)a 3 



2i?x + 45i + -fi 2 S 2 - -P2TR2, (37a) 
P 



P 



(1 - r)a 3 - ^-q(q + 1)(1 + r)a a + U - 1 - i- 9 (g + 1)(1 - T )j a 3 

= 2i? 3 - 45 3 - -/3 2 S 2 + -0 2 tR 2 . (37b) 
p p 

Equations (|37a[) - (|37b)) will be referred to as the reduced Maxwell-like system for & q ^[(j)^]. Given 
a solution 01, a 3 to the reduced system (|37ap - (|37b|l . the coefficient a 2 is obtained from equation 
(|36|) by means of an algebraic manipulation, while ao and 04 are obtained as a solution of the 
ordinary differential equations (|3f a[) and (|3f cp . 



In order to ease the subsequent discussion, the system ()37a[) and (|37b|) is written in matricial 
form as 

y'(T) = A(T)y(T)+b(r), (38) 

with 



Mr) 



1 -p+ ^9(9+ l)(l + r) 
2p I — T 



1 



and 



with 



1 - T 
1 



1 , A-t 

V (<Z + 1) TT7 

p-l-i«(? + l)(l-r) 



y(r) 



01 (r) 

fl 3(T) 



6(r) 



1 + T 



1 



*i(r) 

:F 3 (t) 



i*i = 2iJi + 45i + -P2S2 - -P2TR2, 
P P 

F 3 ee 2R 3 - 45 3 - -P2S2 + -P2TR2. 

P P 



It follows from the discussion of section (|7.2[) that 

F 3 (t)=-F 1 (-t) = -F{(t) 



Following the ideas of the discussion in [53], it is possible to find a fundamental matrix for 
the system (|3"5)) . One obtains: 



Xp, q — 



Qi (-i) q+1 Q 3 
(-i) q+1 QI Qi 



with 



I — r 



P+i 



(hir^(^) P^f'^ir), (39a) 
(^7-)'" * P^r^ir), (39b) 
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where P„ (t), n a non-negative integer, denotes a Jacobi polynomial — see e.g. for defini- 
tions and properties. Furthermore: 

QU^^Q^-t), Ql( T )=Q 3 (- T ). 

The determinant of the fundamental matrix X p . q (the Wronskian) is given by 

AetX lhq = W (1~t 2 )p-\ (40) 

with Wo a constant. Consequently, the inverse X~ q is given by 



/ (1 + r) 2 (-!)« \ 

(I-t)p- 1 M (l-r)f- 1 p '« 



-m . ' , iv P 



where if, L, M and A~ are shorthand for the following Jacobi polynomials 



^ M ( T )^-l T P^+M-rt ( _ T)) 



_ f ,(-p-i,P-i) / 
_ f ,(-p-i,P-i) / 



The solution to system (|38|) is given by 

y(r) =X(t)X- 1 (0) 2 /(0)+X(t) / ^(s^ds. (41) 



o 

One can verify that the procedure described in the previous lines does indeed provide a solution 
to equations (|31a|) - (|31e|) and (|32ap ~ (|32cp . The argument is a follows: 

(i) One solves the Maxwell-like reduced system (|37a|l - (|37bjl using formula (|4ip or any other 
method. 

(ii) Next, one substitutes the values of the coefficients a\ and 03 obtained in this way into the 
evolution (|34b|) . This equation can be solved for a-i by a direct integration. 

(iii) Given the evolution equations (|37a[) . (|34b[) and (|37b|) one can produce an argument to show 
the propagation of the constraint equation ([35]) — if ([35]) is satisfied initially for r = 0, then 
it is satisfied at later times. 

(iv) Equations (|34bj) and ([35| imply the algebraic condition (f36|) . The latter, together with 
(|37ap and (|37bp imply the evolution equations (|34a|) and (|34c|) . 

(v) One substitutes the coefficients 01, 02 and 03 into the Bianchi propagation equations (1 3 lap 
and (|31e|) . Again, these equations can be solved for the coefficients do and 04 by means of 
a direct integration. 

(vi) It can be shown that the evolution equations (|31ap - (|31e|l imply the propagation of the 
constraint equations (|32aP and (|32c|) . 
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7.4 General properties of the solutions to the reduced Maxwell-like 

system 

Let C w (a,6) denote the set of analytic functions on the interval (a, b). From formula |4T|) one 
obtains the following result. 

Proposition 2. If the components of the vector b(r) are polynomials, then the solutions ax, 13 
to the Maxwell-like reduced system 138\) will be either polynomial in r or of the form 

«i = P(r) + ^rQi(r) ln(l - r) + -L(_l)«+iQ 3 ln(l + r), 

«3 = P s (r) + -L(-l)«+iQ-( T ) ln(l - r) + ^Qf(r) ln(l + r). 

where P(t) is a polynomial in t and Q\{t) and Qs(r) are as in H39a\) and \39b\) . These last 
solutions are of class C w (-1, 1) n C p_1 [-1, 1]. 

Proof. The crucial observation is to note that the particular form of the Wronskian ([40]) 
implies that the partial fraction decomposition of the entries in the integrand, X~^b, of formula 
(|4T|) contains negative integer powers of (1 ± r). The terms (1 ± t) _1 will integrate to In |1 ± t|. 
The terms with (1 ± r)~ r , r > 2 integrate to terms of the same form. Multiplication by the 
matrix X Pt g removes these rational terms. □ 

Remark 1. In order to have polynomial-only solutions, some special cancellations should 
occur in the partial fraction decomposition of the entries of X~^b. In the sequel it will be shown 
that these cancellations do occur for particular combinations of the mutiindices (p; q, k). 



Remark 2. The particular form of the Wronskian (|40|) is the reason why the approach 
described in this section has been preferred to the one originally described in [15] . The Wronskian 
of the reduced system advocated in the aforementioned reference is of the form 

C /(t)(1-t)^ 2 , 

with c a constant and 

/(r) = 2(p + l)(p - 1) - (q - l)(q + 2)(1 - r 2 ). 

The partial fraction decomposition of the entries of the corresponding matrix product X~^b will 
contain terms of the form 

ar + f3 



2(p+l)(p-l)-(g-l)(« + 2)(l-T*)' 
for a, (3 some constants. These will integrate to give multiples of terms of the form 



ln|2(p+l)(p-l)-(q-l)(g + 2)(l-T 2 )|, arctan' ' ^ 1 ^ q + 2 ^ 



2(p+l)(p-l)-(?-l)(g + 2) 



As a result of proposition^ one is expecting solutions consisting only of polynomials and In |l±r|. 
Thus, there are some non-trivial cancellations in formula (|4ip that need to be explained by means 
of some further arguments. 



7.5 The Case p = q and the regularity condition at i 

As first pointed out in [15j , for p > 2 if q = p, then the vector b appearing in formula (|4ip is such 
that © Pi fe[6] = 0, so that the solution to the reduced system is given entirely by the solution 
to the homogeneous problem. From here, it is possible to identify conditions on the initial data so 
that the solutions to the Bianchi transport equations for these particular sectors extend smoothly 
through the critical sets X ± . In particular, one has the following result — theorem 8.2 in [15] . 



24 



Theorem 1. Given vacuum initial data which is time symmetric and analytic in a neighbourhood 
of infinity, the solution to the regular finite initial value problem is smooth through only if the 
condition 

D (EpFp ---D ElFl b ABCD) {i)=0, p=l,2,.... (42) 

If this condition is violated at some order p' , then the solution will develop logarithmic singularities 
in & pl [cj)^')] atl ± . 

In the previous result bABCD denotes the spinorial counterpart of the Cotton tensor. The 
Cotton tensor is given by 

kpij = Djl^p D^ljpj lij = Sij -\- -^rhij , 

where s^ denotes the tracefree part of the Ricci tensor of the initial 3-metric . Associated to 
k pi j one has the spinor ItABCDEF, which because of the symmetries of the tensorial counterpart 
can be written as 

kABCDEF — OABCEtDF + OABDF^CE- 

Remark. Time symmetric data which is conformally flat in a neighbourhood of infinity 
satisfies the condition (|42l trivially to all orders. 

8 The solutions to the transport equations for data which 
is Schwarzschildean up to a certain order 

In this section the main analysis of the present work is presented. The key idea behind is to 
analyse the solutions to the transport equation at spatial infinity for initial data sets which are 
Schwarzschildean up to a certain order. 

8.1 Summary of the results of reference |30j 

In order to motivate the analysis of this section, a brief summary of the results of |30j is presented. 
In that reference one assumes a function W of the form 

oo 2p 
p=2 k=0 y ' 

In fact, the form of W assumed in [30J is slightly more general than this as it includes terms 
w i,2,ki fc = 0, 1, 2, which — as seen in section|3] — can always be removed by choosing the centre of 
mass properly. Using scripts in the computer algebra system Maple V one can calculate the ex- 
plicit solutions to the transport equations (|27ap - (|27cp via the decomposition in terms of spherical 
harmonics discussed in section [7] The solutions for the orders p = 0,1,2,3 have been calcu- 
lated in [15j I20j. These solutions have polynomial dependence in r, and thus, extend smoothly 
through the critical sets, although as seen in proposition [2 the solutions could have had loga- 
rithmic singularities. The use of computer algebra methods allows to go beyond this point and 
to calculate further orders in the expansions. For p — 4 one finds again that the solutions have 
polynomial dependence on r. However, for p = 5 the situation is different. One finds that the 
sectors 62 [4>^ 5 '] has solutions with logarithmic terms of the form given by proposition [2j An 
important observation is that these logarithmic solutions do not appear if one chooses a function 
W for which 

W 2 A,k = 0, fc = 0,...,4. 

In the terminology of section [3] this means that the data is Schwarzschildean up to order p = 2. 
Assuming that this is the case, one can proceed further with the expansions. At order p = 6 one 
finds logarithmic solutions only in the sectors 63 [0' 6 ']. The logarithmic solutions can be avoided 
by considering data such that 

w 3 ,6,k = 0, fe = 0, ...,6, 
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that is, data which is Schwarzschildean up to order p = 3. From these results one can already 
infer a pattern which has been confirmed to all the orders for which the calculations have been 
carried out. The calculations reported in [3D] are carried out up to order p = 9, but there is no 
reason why — besides computing power — the calculations cannot be carried out any further. 

The pattern inferred from the analysis in |30j is as follows. Assume that one has an initial data 
set which is Schwarzschildean up to order p = p m . Then, the solutions to the transport equations 
for p < p. will only contain the sectors So[«^]. These will coincide with the sectors implied by 
the s-jet J^'\u^\ of the Schwarzschild spacetime. At order p = p. + 1 the non- vanishing sectors 
are 

Both of them arc polynomial in r. At order p = p. + 2 the non-vanishing sectors are 

6 [ u (p - +2) ], e p . +1 [u^+% & p . +2 [u^'+% 

Again, all the non-vanishing sectors are polynomial in r. At order p = p. + 3 the non-vanishing 
sectors are 

&o[u^+% & p . +1 [u^ + % 6 p . +2 [u^+% 6 p . +3 [u^' + % 

All the non-vanishing sectors are polynomial in r. Finally, at order p = p, + 4 one has the 
following non- vanishing sectors: 

e [u^+% 6 p . +1 [u^+% 6 p . +2 [u^+% & p . +3 [u^+% & p . +4 [u^ + % 

The sectors & [u&' + % 6 p .+ 2 [u (p ' +4) ], 6 p .+ 3 [w (p ' +4) ] and 6 p .+ 4 [u (p,+4) ] have polynomial de- 
pendence in r, while & p ,+i[u( p ' +4 )] will have logarithmic singularities of the type indicated in 
proposition [21 Furthermore, the solutions will be logarithmic free if 

w P .+i,2(p.+i) : fe = 0, k = 0, .. . ,2 (p. + 1). 

This pattern will be effectively proved in the sequel. It readily suggests an inductive procedure 
to prove the main theorem presented in the introductory section. The calculations in [30] are the 
base step of this inductive procedure. 

8.2 Properties of data which is Schwarzschildean up to order p — p. 

We start with some generic observations which will be used systematically in the sequel. Assume 
that the function W , appearing in expression ((4]) for the conformal factor d, is of the form 

oo 2p 

w = j+ E E^wVX' ( 43 ) 

p=p.+i k=a p ' 

so that the initial data is Schwarzschildean up to order p, . The function W — and hence also the 
coefficients w p; 2 Pl /c, p. < p, k = 0, . . . , 2p — appears non-linearly in the expression 



{1 + pWy 

for the conformal factor f2 and, moreover, in the expressions for the initial data for the spinors 
4>abcd and Qabcd on C a;K — see equations (|22a|) and (|22bj) . Hence, when calculating the normal 
expansions of 4>abcd and Qabcd one encounters products of the form T ii 3l ki x T i2 32 k2 with 
■^ii"' I fe ' r ^i2 J2 k2 ^ ^0 w hich have to be linearised — that is, expressed as a linear combination 
of other functions T i J k . This linearisation procedure is extremely cumbersome and involves the 
use of the Clebsch-Gordan coefficients of SU(2,C). Remarkably, it turns out that if one only 
considers expansions up to order p = p. + 4 — which is what will be required in the present 
analysis — these higher order products do not arise. An inspection renders the following result. 
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Lemma 4. For initial data sets which are Schwarzschildean up to order p = p. in B a one has 
that on C a ,K 

4>abcd - Pabcd = 0(P P ' +1 ), Qabcd - O'abcd = 0(p p '+ 2 ), 

where 4>abcd an d ®abcd denote, respectively, the Weyl and Ricci spinors of the Schwarzschild 
data. Moreover, if p, > 3, then 

j 2p.+2 

^ = ^ + (v +1V ^ ™P. + h2p.+2,kT 2p , +2 k pt + 1 _ lP P ' +1 

^* >' k=0 

^ /2p.+2 2p.+4 \ 

T^M I E W P.+^P.+2*T 2pt+2 p >+3 _ j; + E W p . +2 ;2 P .+4,k T 2p.+4 p.+4-i J P P '~ 
>' \ k=0 k=0 I 

2p.+4 2p.+6 \ 

E Wp. + 2:2p.+4,feT' 2p>+4 p>+4 _i + E w P.+3;2p.+6,k T 2p.+e p.+5-i ] P*' 

k=0 k=0 / 



+ 



(p.+ 


2)1 


1 




(p. + 


3)! 


1 




(p. + 


4)! 


o( P p - 


+5) 



/2p.+6 2p.+8 \ 



+4 



where 



Wp^+i-2p,+2i,k — Cp % +i-2p.+2i,kU!p t +i-2p,+2i,k, 

with c p .+i ; 2p.+2i,fc some numerical coefficients which can be explicitly calculated, and 

4>\ = —6m, i = 2, 
« = 0, *^2. 

j4 similar expansion holds for the components of Oabcd ■ 

The following result, which can be proved by direct inspection, shows that the calculations 
described in the present work do not require the calculation of complicated SU(2, C) Clebsch- 
Gordan coefficients. 

Lemma 5. For the class of initial data under consideration, the terms 

E( p ) (q(»w,«^) + lC 



3=1 

p 



-.O'UCp-j) 



J'=l w 

E ( P ) (BF^KM-F^idXjMa^) 



3=1 

with p = p, + 1, . . . ,p, + 4 cmc? p. > 3, appearing, respectively, in the transport equations \27aty , 



(27b) and \21c^ contain, before linearisation only products of the form T x T { \ 



A direct observation is the following. 

Lemma 6. Let [v] and [4>] be the s-jets of order p arising from initial data with ADM 
mass m which is Schwarzschildean up to order p = p % . Then one has that 

for p = 0, . ..,£>,, where J% [v m ] and Jj P '[</>•] are the s-jets of order p arising from exactly 
Schwarzschildean data with ADM mass m. The first difference between these two sets of jets 
arises at order p = p % + 1 . 
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Remark. An inspection of the terms discussed in lemma [5J for p = p, + 4 shows that the 
present discussion requires at most the explicit knowledge of the s-jet [u,] of the Schwarzschild 
spacetimc. 

In order to appreciate the following arguments, the non- vanishing sectors in the s-jets J% +5 ^ [v- 
v,] and j!f' +jS> [4> — 4>m\ w iU be listed. This list can be deduced from lemma[TJ 

• At order p = p, + 1 : 



At order p = p, + 2: 



• At order p = p, + 3: 

6 [v<**+% 6 p . + M p ' + % 6 P .+2[v (p ' + % 

6 [^ p ' + % e p . +1 [^-+% & p . +2 [^-+% 6 p . +3 [^' +3) ]. 

• At order p = p, + 4: 

6 {v (p ' + % e p . +1 [v^+% 6 p . +2 [ v ( p ' + % 6 p . +3 k> +4) ] 
&ott (p ' + % e p . +1 [^- + % & p . +2 [^- + % & Pm+s [^' + % s p . +4 [ 



(P. +4)1 



Remark. From lemma ([5]) it follows that there is no mixing between the various sectors arising 
at each order. Thus, it is only necessary to carry out a discussion of the solutions in the difference 
jet Jj P * +5 ' > [u — u,] for the sectors 6p.+i. The analysis of the corresponding sectors 6p.+2, & P .+3, 
S p .+4 and & P ,+5 can, in principle, be obtained from that of the sector ©p.+i by performing, 
respectively, the formal replacements p, i— *p. + 1, p% i— > p, + 2, p m > p. + 3 and p. i— > p. + 4. 

Warning. In order to improve the readability, obvious strings of subindices will be omitted in 
the sequel. 

8.3 Properties of the transport equations for p — p. + 1 

The solutions to the transport equations at order p = p m + l can be essentially read from the orig- 
inal analysis carried out in [15] . Using lemma ([I]) it follows that t/P« +1 ) contains no contribution 
to the sectors 6 p .+i. The only non-vanishing sector in is So, which coincides with the 

Schwarzschildean solution. The (non-vanishing) contribution of to the sectors & P .+i is 

given by: 

ao )P .+i;2(p.+i),* = A),fc(l ~ t)p-+ 3 (1 + t)p--\ (44a) 
ai, P .+i;2( P .+i), fe = A 1)fc (l - r)P'+ 2 (l + rf , (44b) 
a 2 , P . + i;2(p.+i), fc = A 2 , k (l - rf' +1 {l + rf'+\ (44c) 
a 3 , P .+i;2( P .+i),k = Ai, fc (l - r) p - (1 + r) p '+ 2 , (44d) 
a4, P .+i;2( P . + i), fc = A),fc(l ~ rf-^l + rf'+\ (44e) 
with A: = 0, . . . , 2 (p. + 1) and 

A),fc = -y/P*(P» + IXP* + 2 )(P» + 3 ) w p. + l;2(p. + l),fc, 

Ai ifc ee -4(p. + 3)\/ (p. + l)(p. + 2)w p . + i ;2 (p.+i),fe, 
^2,fe = -6 (p. + 2) (p. + 3)w p . +1 . 2 (p. + i), fc . 

Hence, the solutions extend analytically through the sets X^, consistently with the theorem 
concerning the regularity condition (|42p . 
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8.4 Properties of the transport equations for p = p, + 2 

The solutions for & p ,+i[v^ p ' +2 ' ) ] can be calculated using the solutions for & p .+i [</>' p,+1 '] given by 
(I44a[) - (|44cp . According to lemma[2]the solutions for & p ,+i [i/ p,+2 - ) ] will be polynomial in r and it 
turns out that they can be expressed in terms of hypergeometric functions. With this information 
at hand it would be, in principle, possible to analyse the Maxwell-like reduced system (|37aj) - (|37bj) 
order p = p, + 2 for the coefficients ai. p .+2 and a3 iP# +2- However, the form of the explicit solutions 
in 6 p#+ i[i/ p,+2 )] make it very difficult to identify useful structures in the equations. Instead, it 
is desirable to find a way to extract information on the solutions to the p — p m + 2 MaxwelHikc 
reduced system without having to make use of the explicit solutions of 6p >+ i[t/ p * +2 )]. 

In order to get around this problem we do the following: a close inspection of the right hand 
side of equations (|37a|) and (|37bj) reveals that by differentiating three times with respect to r and 
then using the u(p»+ 2 ) and transport equations and r-derivatives thereof, it is possible to 

obtain a reduced system for 

[3] _ ^3 [3] _ o3 

a l,p.+2 = °r a l,Z>.+2> a 3,p.+2 = °V a 3,p.+2, 

where the non-homogeneous terms depend explicitly only on ai^.+i and 03^+1 as given by 
expressions (|44b|) and (|44dj) . Clearly, if the solutions a^ Pt+2 and a 3 3 | )>+2 for these equations are 

regular at r = ±1, this will also be the case for ai iP>+ 2 and 03^+2- Similarly, if Qj p <+2 and 

[31 

a 3p.+2 have logarithmic singularities at r = ±1, also ai. p>+2 and d3 iP> +2 will also contain them. 

[31 [31 

The calculation of suitable equations for a[ Pt+2 an< i a 3 p. +2 requires systematic and extensive 
use of computer algebra methods. Calculations in the system Maple V render a reduced system 
of the form: 

(1 + r)d T a l ^ +2 + ({p. + l)(p. + 2)r - {p\ - 9p. + 2))ag,. +2 



1 



2( P .-1) (P ' + + 2)(1 ~ T ) a 3 3 L + 2 



rGp. +2 ai, p .+i + j- .„ . — —H p , +2 a 3 ^, + i 7 (45a) 



(l-r) 3 (l + r) 4 (l-r) 3 (l + r) 4 

(1 - r)d T af pt+2 - ^^ (P. + 1)(P. + 2)(1 + r)a [ ^ +2 



1 



2(P. - 1) 



(p. + l)(p. + 2)r+(p 2 -9p. + 2 
1 



J3] 



- (1 + T )3(1 _7)4 i7 p.+2 a l.P. + l + (JTP7)37J _ T )4 G P.+2 a 3,P. + l- (45b) 

In the above expressions Gp.+2 and H p , +2 are explicit polynomials in r of degree 9 with coeffi- 
cients which are themselves polynomials in p.. Furthermore, G pt+2 (r) = Gp.+2(~ t)) H p ^ +2 (r) = 
H Pm +2(—T~). Both G p .+2 and i?p.+2 contain an overall factor of m. The explicit form of these 
polynomials is given in the appendix. The system (I45aj) - (|45bj) is supplemented by the initial 
conditions 

4 3 L+2(0) = -Wi(6p. + 23)(p. + 2), 4 3 L+2(°) = 2mA! (6p. + 23)(p. + 2). (46) 

The data is calculated by using the values of <^ p ' +1 )(0), (p ' +2) (O), z/ p>+2) (0) implied by the data 
jets Jt'" +2 \v], J^' +2 \4>\ and by r-differentiating as necessary the u(p«+ 2 ) and ^(p«+ 2 ) transport 
equations. This cumbersome calculation has been carried out in the computer algebra system 
Maple V. 

Substitution of the explicit tfi( p,+1 ' > solutions (|44a[) - (|44e[) into equations ()45a[) and (|45b[) shows 
that the right hand sides of these equations are, respectively, of the form 

(1 - t)*"" 3 (1 + Tf-- 4 Qp. + 2(T), (1 + r)f- 3 (l - rf -^Q^^t), 
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where Q p ,+2{t) is a polynomial of degree 11 in r such that Q p , +2 (±l) 7^ 0. Consistently with 
the above, it is assumed that p, > 4. The cases with p, < 4 can be analysed in a case by case 
basis — cfr. the calculations in 30J . 

The remarkable structure of the zeros of the right hand sides of equations (|45a|) and (|45bj) 
eases the task of looking for polynomial solutions to these equations. Indeed, we note the following 
lemma. 

Lemma 7. All polynomial solutions of equations |^5a[ ) and \45b\) for p m > 4 are of the form 

41 +2 (t) = (l-r)"-- 2 (l + rf-\ +2 (r), (47a) 
41 +2 W - -(1 + - rY'-%l +2 (r), (47b) 

where & p .+2(t) is polynomial of degree 9 and b pt+2 {T) = & p ._|-2(— t). 

Proof. The proof of this lemma is inspired by the discussion in chapter 4 of [23) . One starts by 
looking at the possible zeros of the solution to equations (|45aP and (|45bp at r = 1. If aij 3 ' and 

[31 

a!; are polynomial then one can write 

n* m* 

of =J] ttt (l- r) fc , 4 31 - £ " r) fe , (48) 

k— n* k—m^ 

for some integers n*, m*, n*, m* > and some complex numbers ak, Pi, n* < k < n*, 
to, < I < m* . Dividing equation ([45ap by 1 + r and equation (|45bll by 1 — r and then substituting 
the expressions (|4"5|) into the ordinary differential equations ([45ap and (|45b[) one finds that 

mm{n* — 1, ?7i* + 1} = p, — 3, min{n* — 1, to, — 1} = p. — 5, 

as r = 1 is a zero of the right hand sides of (|45a[) and (|45bp and, moreover, the left and right 
hand sides must have the same multiplicity. The above conditions are satisfied by setting 

n* = p, — 2, to, = p. — 4. 

The discussion of zeros at r = — 1 follows by symmetry — cfr. the discussion in section l?T2"l The 
degree of the polynomial 6 p .+2(t) follows by inspection of Q p ,+2- d 

Now, the substitution of the Ansatz (|47ap and (|47bp with 

9 

fc=0 

into the equations ()45a|) and (|45b[) leads to a system of 11 linear algebraic equations for the 10 
unknowns i? p .+2,fc: k — 0, ... ,9. This overdetermined system can be seen to have a solution which 
can be explicitly calculated. The result of this calculation is also presented in the appendix. In 
particular, the solution obtained is such that 

«i 3 L+ 2 (0) = -2mAi(6p. + 23) (p. + 2), ajg, #+2 (0) = 2mA 1 (6p. + 23) (p. + 2), 

consistent with the initial data (|4T)]) for the system (|47a|l and (I47b|) . Thus, the particular solution 
to the system (|45al) - (|45bp calculated by the above procedure is, in fact, the solution to the system 
in question with data given by (|46p . 

From the procedure described in the above paragraphs together with explicit calculations in 
[50] for p. < 4 one has the following result. 



Proposition 3. The solution of the Maxwell-like reduced system |^5a[ ) and (J^5b) with data given 
by |^g[ ) is polynomial in r and thus it extends smoothly through r = ±1. 

[31 [31 

Remark. It follows that not only a[ Pt+2 anc ^ a 3,p.+2 but also &i.p >+ 2 an d Q3,p.+2 are polynomial 



in r. Moreover, from the discussion in section [7. 31 and lemma[7]also ao iP .+2, a2,p.+2, fl4.p.+2 and 
the sector & p .+i [u' p,+3 - ) ] are polynomial in r. 
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8.5 Properties of the transport equations for p = p, + 3 

We extend the approach used in the previous section to decide whether the coefficients ai !P# +3 
and d3 iP> +3 are polynomial in r or not. For these coefficients, the calculations are more involved, 
and require r-differentiating the Maxwell-like reduced equations satisfied by ai jP .+3 and 03^+3 
seven times to obtain a linear system for 

[7] _ A 7 [7] _ A 7 

a l,p.+3 = °T a l,Z>.+3> a 3,p.+3 = °V a 3,p.+3- 

After computer algebra computations involving the substitution of the transport equations satis- 
fied by t>(P«+ 2 ), y(P'+ 3 ) and cf>( p,+1 \ 0(p»+ 2 ) and r-derivatives thereof, one obtains the following 
system of equations: 

(1 + r)d T a [ ^ +3 - — L ((p. + l)(p. + 2)r - (p 2 . - 21p. - 38))ag,. +3 

~ 2(p.-4) (p> + 1)iP ' + 2){1 ~ T)a S"+ 3 

1 _ 1 „ 

- (l-r)8(l+r)9 Gp - +3ai ^ +1 + (l-r)8(l + r)9 Hp - +3fl3 ' ?, - +1 

1 1 K [3] 1 j [3] 

+ (1 - T) 4 (l + T )5 A P-+ 3a l^.+2 + (1 _ r )4 (1 + T) 5^P.+3 a 3,P.+2 J 

(49a) 

(1 - r)d T ag,. + 3 + 4) (P. + 1)(P. + 2)(1 + r>' 7 p . + 3 

> ( >.^4) { {p ' + 1){p ' + 2)r + {p2 ' ~ 2lPm + ^^H"-*- 3 



(l + r) 8 (l-r) 9 (1+t) 8 (1-t) £ 



1 T s n [3] 1 1 V s r, [3] 



T (l + T) 4 (l - r )5^P.+3"l,P.+2 T (1 + t) 4 (1 _ r) 5-p.+3"3, P .+2- 

(49b) 

In the above expressions G p> +3 and iJ p>+ 3 are explicit polynomials in r of degree 19, while 
Kp.+3 an d ^p.+3 ar e explicit polynomials of degree 10. The polynomials G p ._|_3 and -Hp. +3 contain 
an overall factor of m 2 while K Pt+ 3 and £ p .+3 have an overall factor of m. The coefficients of 
these polynomials are themselves polynomials in p.. As in the case of p = p. + 2, initial conditions 
for the system (|4^al)- p9b)) are calculated using the values of ^P' +1 ~>(0), <^' +2 )(0), ^P'+ 3 ) (0), 
w (p.+2)( ), w (p.+3)(o) implied by the data jets J^* +2) M, Jj'^M- 4 P ' + %L 4 P * + %] and 
by r-differentiating as necessary the ?j(p»+ 3 ) and 0(p«+ 3 ) transport equations. One obtains 



< P .+ 7 (0) = 77™ 2 ^(p. + 2)(p. + 3)(3p5 + 7 2p4 + 49353p 3 



8 

+ 142260^ + 610272p. + 302048), (50a) 
4 7 L+7(0) = -\™ 2 A{ P . + 2)(p. + 3)(3^ + 7 2p4 + 49353p 3 

„2 



+142260^ + 610272p. + 302048). (50b) 



We proceed to analyse this system on the same lines as it was done for (|45a|) - (|45b[) . First one 
notices that direct substitution of the explicit solutions ai )P .+i, c^p.+i, ai iP .+2 and 03 iP .+2 which 
were obtained, respectively, in sections 1831 and I8T41 shows that the right hand sides of equations 
(|49a|l and (|49b|) are of the form 

(1 - tT-- 8 (1 + rr- 9 Q p . +3 (r), (1 + rf- 8 (l - r)^- 9 Q; +3 (r) 
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where Q p ,+s(t) is a polynomial in r such that Q p .+3(±1) =/= 0. In order for the following 
calculations to make sense, it is assumed that p, > 9. The cases p. < 9 can be analysed 
individually and provide the same qualitative picture. 

The remarkable structure of the zeros of the right hand sides of equations (|45ap and (|45b[) 
eases the task of looking for polynomial solutions to these equations. Indeed, one has the following 
lemma, whose proof is similar to that of lemma [Jj 

Lemma 8. All polynomial solutions of equations \49a\) and (49b\ ) are of the form 

«S. +3 (r) = (1 - rf- 7 (l + Tf-%. +3 (T), (51a) 
41 +3 W = -(1 + t) p '- 7 (1 rr- 9 bl +3 (r), (51b) 
where 6 p .+3(r) is polynomial of degree 19, and bp t+3 (r) = 6 p .+3(— r). 

Now, we use the expressions (|51a[) and (|51bj) as an Ansatz for oJP Pt+ 3 and a3 7 p>+3 with 

19 

The substitution of this Ansatz into equations (|49a|) - (|49bj) leads to a system of 21 linear algebraic 
equations for the 20 unknowns B Pt+ 3^, k — 0, . . . , 19. By means of an explicit calculation, this 
overdetermined system can be seen to have a solution — that is, not all the 21 equations are 
linearly independent. The solution so obtained contains a free parameter which can be adjusted 
to set the value of the coefficient B p ^ + ^ : q such that the initial conditions (|50a|) - (|50bjl are satisfied. 

With the procedure described in the above paragraphs and explicit calculations for the cases 
with p, < 9 one has the following result. 

Proposition 4. The solution of the Maxwell-like reduced system \49aty - l49ty with data given by 
i50a\) - [50b^) is polynomial in r and thus, it extends smoothly through r = ±1. 

[71 [71 

Remark. It follows that not only a[ Pt+ 3 and (%L +3 but also ai iP .+3 and a3 jP .+3 are polynomial 
in r. Moreover, from the discussion in section [7731 and lemma ao. p .+3 a2, p .+3, a4, P .+3 and the 
sector 6p.+i [i/ p * +4 - ) ] are polynomial in r. 

8.6 Properties of the transport equations for p = p, + 4 

Finally, it is shown, by methods similar to the ones used for the cases p = p, + 2 and p = p, + 3 
that the solutions ai iP> +4 and 03^+4 of the order p = p % + 4 MaxwelUike transport equations 
cannot be purely polynomial. In order to show this, one r-differentiates eleven times(!) the 
MaxwelUike reduced equations satisfied by ai jP .+4 and 03^.4.4 to obtain a linear system for 

[11] _ oil [11] _ all 

a l,p.+4 = °t a l,p.+4, a 3, p .+4 = °t a 3,p.+4- 

After lengthy computer algebra computations involving the substitution of the transport equa- 
tions satisfied by v p,+4 , v p ' +3 , v p,+2 , P,+3 , </> p,+2 , and r-derivatives thereof into the 
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equations for 4 1 pl+4 and 4p.-i-4' one obtains a system of equations of the form: 

[11] 
p. 

1 .. . _ J- 1 J 



(1 + r)d T a[^ +i - _ [{p. + l)(p. + 2)r - (pi - 33p. + 110) )a l $ m+i 



2(p> _ 7) (P. + 1)(P. + 2)(1 - r)a^ J . +4 

= (1 - r) 13 (l + T )i4 G P.+ 4ai ^.+ 1 + (1 - T )i3(i _|_ r )i4- ff P.+ 4a 3,P.+i 

! ^ [3] 1 r [3] 

-Ap. +4 a lif)>+2 + - — r^— ^p. +4 a 3!p>+2 



(l-r)°(l + r) 10 "•"^ (l-r) 9 (l+r) ] 

1 [7] 1 [7] 

+ (l-r)4(l + r)5 Mp - +4ai ^+ 3 + (l-r)4(l + r)5 Afp - +4a3 ^+ 3 ' 

(52a) 

(1 - r)d T a [ ^ +4 + _L_-(p. + l)(p. + 2)(1 + r) a M +4 

~ 2(p.-7) ( (p ' + + 2)T + (P * " 33P * + 11 °)) a £ I .+4 

= (1 + r) 13 (l - T )14- ff P.+4 ai >P' + 1 + (J + T )13(l _ T )14 G P.+4 a 3,P. + l 



(1 + t) 9 (1 - r )10 P.+ 4 "l,P.+2 (1 + r )9(l _ T )10' fl P.+4 u '3,p.+2 

+ (1 + T)4(l - r)5 J V+4 a l,p.+3 + (1 + t)4(1 _ r) 5 M P.+* a 3,p.+3- 

(52b) 

In the above expressions G p> + 4 and H Pm+ 4 are polynomials in t of degree 29 with an overall factor 
of m 3 ; iip >+4 and £ p>+4 are polynomials in r of degree 20 with an overall factor of to 2 ; finally, 
M p>+4 and A r p>+4 are polynomials in t of degree 9 with an overall factor of m. The coefficients 
of all these polynomials are themselves polynomial in p,. The initial conditions for the system 
(j52a|) - (l52b|) can be calculated following a similar procedure to the one used for the lower order 
systems (I45a|) - (|45b|) and (I49a|) - (|49b|) . The initial conditions read 

a i 1 p 1 .+ 4 ( ) = -^m 3 A(p. + 2)(p. + 3) (p. + 4)(1737p£ + 64317^ + 84931161^ - 175033869p 4 
+3959957166p 3 - 1379005296p 2 + 6565457856p. + 1464602048) (53) 

4*1 +4 (°) = l m3A (P- + 2) (p. + 3) (p. + 4)(1737p^ + 64317p^ + 84931161p^ - 175033869p 4 
+3959957166p 3 - 1379005296p 2 + 6565457856p. + 1464602048). (54) 



As in the previous cases, the crucial observation is that using the explicit expressions for 

ai,p.+i, a 3,p.+i> a i 3 |,.+27 4 3 p.+2' a i 7 p.+3 anc ^ a 3p.+3 discussed in the previous sections, one finds 
that the right hand sides of equations (|52ap and (|52b[) are, respectively, of the form 

(1 - rf- 13 (l + Tf- 14 Q p . +4 (r), -(1 + - 13 (1 - r^- 14 Q^ +4 (r), 

with Q p ,+a{t) a polynomial in t such that Q Pe _|_ 4 (±l) 7^ 0. In what follows it is assumed that 
p. > 14, so that the discussions make sense. The cases p. < 14 can be dealt with direct case-by- 
case calculations. The qualitative results are the same as in the general case. 

As in the cases discussed in sections 18.41 and 18.51 the structure of the polynomial solutions (if 
any) of the system (|52a[) - (|52b[) is very particular. More precisely, one has the following result. 

Lemma 9. The polynomial solutions (if any) of the reduced system h52a\) and \52b\) are of the 
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form 



"S+iW = (l-r) p - - 12 (l + rf -%. +4 (t), (55a) 
4!?. +4 (r) = -(1 + rr- 12 (l - r)^ 14 5; +4 (r), (55b) 
where fr p .+4 is a polynomial of degree 29. 



As in the previous orders, we use expressions (|55a[) - (|55b[) as an Ansatz for a[J t+i and a^ m+i 
with 

29 

Again, this Ansatz leads to a system of 31 linear algebraic equations for the 30 unknowns -B p .+4 l fe, 
fe = 0, . . . , 29. By means of an explicit calculation, this overdetermined system can be seen to 
have no solutions unless io p .+i,2p.+2,fc = or m = 0. This is the crucial result of our analysis. 
One has the following proposition. 

Proposition 5. The solution of the Maxwell-like reduced system 1152 a)) and H52b\) admits no 
polynomial solutions unless UJp.+i 2p.+2,k = or m = 0. 

This last proposition, together with proposition^ in section 17741 renders the following corollary. 

Corollary 1. The solution of the Maxwell-like reduced system (5Ba)) - f52b)) with data given by 
h53\) -\54ty develops logarithmic singularities at r — ±1 unless u?p.+i,2p,+2,fc = or m = 0. The 
logarithmic solutions are of class C u (— 1, 1) f] C p ' +3 [— 1, 1]. 

Remark. From the discussion in section 17731 it follows that the coefficients ai }Pm +4, i = 0, .. .,4 
will have logarithmic singularities of the form given by proposition [2] 



9 The main result 

The discussion in the previous sections is summarised in the following result. 

Proposition 6. Given a time symmetric initial data set which in a neighbourhood B a of infinity 
is Schwarzschildean up to order p — p., the solutions to the transport equations for the orders 
p = p, + 1, p % + 2, p % + 3 are polynomial in t and hence extend smoothly through the critical sets 
I . On the other hand, the solutions at order p = p, + 4 contain logarithmic singularities which 
can be avoided if and only if the initial data is, in fact, Schwarzschildean to order p = p, + 1 . 

From this result, an induction argument which uses the explicit calculations performed in |30j 
as a base step renders our main result — cfr. the main theorem in the introduction. 

Theorem 2. The solution to the regular finite initial value problem at spatial infinity for initial 
data which is time symmetric and conformally flat in a neighbourhood of infinity is smooth through 
I if and only if the restriction of the data to T° coincides with the restriction of Schwarzschildean 
data at every order. Furthermore, the analyticity of the data implies that the initial data is exactly 
Schwarzschildean in a neighbourhood of infinity. 

As mentioned in the introduction, the evidence gathered in [29j suggests that it is possible 
to obtain a generalisation of this result for initial data sets which are not conformally flat. In 
that case, the expected conclusion is that smoothness at the critical sets implies staticity at every 
order. In [19j it has been shown that the conformal structure of static spacetimes is as regular as 
one would expect it to be. In particular, static spacetimes extend smoothly (in fact analytically) 
at the critical sets X ± . Key to such a generalisation of theorem [2] is to consider a class of time 
symmetric initial data sets for which it is simple to decide whether its development will be a 
static spacetime or not. 
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A Some of the polynomials in section [8] 

In order to exemplify the type of expressions one has to deal with in the analysis described in 
section [51 here we present some of the simpler polynomials. 

A.l The polynomials in section 18.41 

, 1 , G P .+2 = (p. + 4) (p. + 2) (2p. 3 + 10p. 2 - 2p. - 55) r 9 - (l7p. 3 + 635 + 305 p. + 14p. 4 + 2p. 5 + 8p. 2 ) r 8 
m(p. + 2) 

-2 (175 p. 2 + 7p. 5 - 550p. + 246p. 3 - 652 + 72 p. 4 ) r - 2 (5p. 5 + 155p. 2 + 115p. 3 + 34p. 4 - 471 p. - 1080) r 6 + 
(833p. 3 - 1124 + 795p. 2 + 12p. 5 - 1172p. +206p. 4 ) r 5 + (36p. 5 - 2494 + 601 p. 3 - 946p. + 166 p. 4 + 891 p. 2 ) r 4 
+2 (263 p. + 70 + 12 p. 4 - 280p. 2 - 263 p. 3 ) t 3 - 12 (-32 p. - 95 + 78 p. 2 + 12 p. 4 + 22 p. 3 ) t 2 
-3 (-62 p. 2 - 72 + 51p. 3 +44 p.) t - 75 + 45 p. + 54 p. 3 + 39p. 2 , 

' H p . +2 = -2 (p. +4) (p. + 2) (2p. 3 + 10 p. 2 - 2p. - 55) r 9 +2 (161 p. + 102 p. 2 + 18p. 4 + 33 + 53 p. 3 + 2p. 5 ) r 8 



m(p. + 2) 

+4 (p. + 2) (7p. 4 + 56p. 3 + 122 p. 2 - 71 p. - 450) r 7 + 4 (3 p. 3 - 113p. 2 - 305 p. + 20p. 4 - 78 + 5p. 5 ) r B 

-2 (-1352p. + 905p. 3 + 959p. 2 - 2748 + 12p. 5 + 226p. 4 ) r 5 - 2 (245p. 3 + 36p. 5 + 142 p. 4 - 866p. - 246 - 13p. 2 ) r 

+4 (335p. 3 + 24 p. 4 - 359 p. + 562p. 2 - 886) r 3 + 48 (lip. 3 - 7p. - 13 + 6p. 4 ) t 2 

-6 (-176 + 98 p. 2 + 21 p. 3 + 60p.) r + 90 + 78p. - 108p. 3 - 186p. 2 , 

-bp+2 (r) = -2 (6 p. + 23) (p. + 2) - 2 (p. + 2) (6 p. 2 + 14 p. + 5) r + 4 (p. + 2) (21 p. 2 + 45 p. + 73) r 2 



+4 (p. + 2) (7p. 3 + 16p. 2 + 88 p. + 48) r 3 - 2 (p. + 2) (20p. 3 + 55p. 2 + 229 p. + 256) r 4 

-2 (p. + 1) (p. + 2) (4p. 3 + 15p. 2 + 122p. + 214) r 5 + - (p. + 2) (7 p. 3 + 69 p. 2 + 260 p. + 249) r 6 

+ ^ (p. 3 + 29p. 2 + 127 p. + 120) (p. + 2) 2 r 7 - | (p. + 2) (2p. 3 + 27p. 2 + 109p. + 99) r 8 

~ (p. + 2) (2 p. 4 + 23 p. 3 + 97 p. 2 + 178 p. + 111) t° . 

The expressions of the polynomials appearing in the analyses at order p = p. + 3 and p — p, +4 
are too unwieldly to be presented here. 

References 

[1] P. Anninos, D. Hobill, E. Seidel, L. Smarr, & W. M. Suen, Collision of two black holes, 
Phys. Rev. Lett. 71, 2851 (1993). 

[2] A. Ashtekar, Lectures on non-perturbative canonical gravity, World Scientific, 1991. 

[3] R. Beig, Conformal properties of static spacetimes, Class. Quantum Grav. 8, 263 (1991). 

[4] D. R. Brill & R. W. Lindquist, Interaction energy in geometrostatics, Phys. Rev. 131, 471 
(1963). 

[5] P. T. Chrusciel & E. Delay, Existence of non-trivial, vacuum, asymptotically simple space- 
times, Class. Quantum Grav. 19, L71 (2002). 



35 



[6] S. Dain & H. Friedrich, Asymptotically flat initial data with prescribed regularity at infinity, 
Comm. Math. Phys. 222, 569 (2001). 

[7] J. Frauendiener, Numerical treatment of the hyperboloidal initial value problem for the vac- 
uum Einstein equations. I. The conformal field equations, Phys. Rev. D 58, 064002 (1998). 

[8] H. Friedrich, The asymptotic characteristic initial value problem for Einstein's vacuum 
field equations as an initial value problem for a first-order quasilinear symmetric hyperbolic 
system, Proc. Roy. Soc. Lond. A 378, 401 (1981). 

[9] H. Friedrich, On the regular and the asymptotic characteristic initial value problem for 
Einstein's vacuum field equations, Proc. Roy. Soc. Lond. A 375, 169 (1981). 

[10] H. Friedrich, On the existence of analytic null asymptotically fiat solutions of Einstein's 
vacuum field equations, Proc. Roy. Soc. Lond. A 381, 361 (1982). 

[11] H. Friedrich, On purely radiative space-times, Comm. Math. Phys. 103, 35 (1986). 

[12] H. Friedrich, On the existence of n-geodesically complete or future complete solutions of 
Einstein's field equations with smooth asymptotic structure, Comm. Math. Phys. 107, 587 
(1986). 

[13] H. Friedrich, On static and radiative space-times, Comm. Math. Phys. 119, 51 (1988). 

[14] H. Friedrich, Einstein equations and conformal structure: existence of anti-de Sitter-type 
space-times, J. Geom. Phys. 17, 125 (1995). 

[15] H. Friedrich, Gravitational fields near space-like and null infinity, J. Geom. Phys. 24, 83 
(1998). 

[16] H. Friedrich, Conformal Einstein evolution, in The conformal structure of spacetime: Ge- 
ometry, Analysis, Numerics, edited by J. Frauendiener & H. Friedrich, Lecture Notes in 
Physics, page 1, Springer, 2002. 

[17] H. Friedrich, Conformal geodesies on vacuum spacetimes, Comm. Math. Phys. 235, 513 
(2003). 

[18] H. Friedrich, Spin-2 fields on Minkowski space near space-like and null infinity, Class. 
Quantum Grav. 20, 101 (2003). 

[19] H. Friedrich, Smoothness at null infinity and the structure of initial data, in 50 years of the 
Cauchy problem in general relativity, edited by P. T. Chrusciel & H. Friedrich, Birkhausser, 
2004. 

[20] H. Friedrich & J. Kannar, Bondi-type systems near space-like infinity and the calculation of 
the NP-constants, J. Math. Phys. 41, 2195 (2000). 

[21] C. Liibbe & J. A. Valiente Kroon, On de Sitter-like and Minkow ski-like spacetimes, 2009. 

[22] C. W. Misner, The method of images in geometrodynamics, Ann. Phys. 24, 102 (1963). 

[23] M.van der Put & M. Singer, Galois Theory of Linear Differential Equations, Springer, 2003. 

[24] R. Penrose, Asymptotic properties of fields and space-times, Phys. Rev. Lett. 10, 66 (1963). 

[25] R. Penrose & W. Rindlcr, Spinors and space-time. Volume 2. Spinor and twistor methods 
in space-time geometry, Cambridge University Press, 1986. 

[26] P. Sommers, Space spinors, J. Math. Phys. 21, 2567 (1980). 

[27] J. Stewart, Advanced general relativity, Cambridge University Press, 1991. 

[28] G. Szego, Orthogonal polynomials, volume 23 of AMS Colloq. Pub., AMS, 1978. 



36 



[29] J. A. Valiente Kroon, Does asymptotic simplicity allow for radiation near spatial infinity?, 
Comm. Math. Phys. 251 (2004). 

[30] J. A. Valiente Kroon, A new class of obstructions to the smoothness of null infinity, Comm. 
Math. Phys. 244, 133 (2004). 

[31] J. A. Valiente Kroon, Time asymmetric spacetimes near null and spatial infinity. I. Expan- 
sions of developments of conformally flat data, Class. Quantum Grav. 23, 5457 (2004). 

[32] J. A. Valiente Kroon, Time asymmetric spacetimes near null and spatial infinity. II. Ex- 
pansions of developments of initial data sets with non-smooth conformal metrics, Class. 
Quantum Grav. 22, 1683 (2005). 

[33] J. A. Valiente Kroon, On smoothness- asymmetric null infinities, Class. Quantum Grav. 23, 
3593 (2006). 

[34] J. A. Valiente Kroon, The Maxwell field on the Schwarzschild spacetime: behaviour near 
spatial infinity, Proc. Roy. Soc. Lond. A 463, 2609 (2007). 

[35] J. A. Valiente Kroon, Estimates for the Maxwell field near the spatial and null infinity of 
the Schwarzschild spacetime, 2009. 



37 



